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Abstract. This is a complement to my previous article “Advanced Determinant 
Calculus” {Seminaire Lotharingien Combin. 42 (1999), Article B42q, 67 pp.). In 
the present article, I share with the reader my experience of applying the methods 
described in the previous article in order to solve a particular problem from number 
theory (G. Almkvist, J. Petersson and the author, Experiment. Math. 12 (2003), 441- 
456). Moreover, I add a list of determinant evaluations which I consider as interesting, 
which have been found since the appearance of the previous article, or which I failed 
to mention there, including several conjectures and open problems. 


1. Introduction 

In the previous article unni, I described several methods to evaluate determinants, 
and I provided a long list of known determinant evaluations. The present article is 
meant as a complement to uni. Its purpose is three-fold: hrst, I want to shed light 
on the problem of evaluating determinants from a slightly different angle, by sharing 
with the reader my experience of applying the methods from uni in order to solve a 
particular problem from number theory (see Sections El and E}; second, I shall address 
the question why it is apparently in the hrst case combinatorialists (such as myself) 
who are so interested in determinant evaluations and get so easily excited about them 
(see Section H; and, hnally third, I add a list of determinant evaluations, which I 
consider as interesting, which have been found since the appearance of |l()9j , or which I 
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failed to mention in the list given in Section 3 of uni (see Section li), including several 
conjectures and open problems. 


2. Enumerative combinatorics, nice formulae, and determinants 


Why are combinatorialists so fascinated by determinant evaluations? 

A simplistic answer to this question goes as follows. Clearly, binomial coefficients 
(^) or Stirling numbers (of the second kind) S{n, k) are basic objects in (enumerative) 
combinatorics; after all they count the subsets of cardinality fc of a set with n elements, 
respectively the ways of partitioning such a set of n elements into k pairwise disjoint 
non-empty subsets. Thus, if one sees an identity such as^ 


dt ( ( ^ ^ ~ ~ !)• 

V V« -^+j) ) (« + * - 1)! (^ + ^ - 1)!’ 


( 2 . 1 ) 


det (^(i+j,i))=JJf (2.2) 

l<2,7<n 

“ i=l 

(and there are many more of that kind; see jlOh] and Sectional), there is an obvious 
excitement that one cannot escape. 

Although this is indeed an explanation which applies in many cases, there is also an 
answer on a more substantial level, which brings us to the reason why / like (and need) 
determinant evaluations. 

The favourite question for an enumerative combinatorialist (such as myself) is 


How many (...) are there? 

Here, (...) can be permutations with certain properties, certain partitions, certain 
paths, certain trees, etc. The favourite theorem then is: 


Theorem 1. The number of {...) of size n is equal to 

NICE{n). 


I have already explained the meaning of (...). What does NICE{n) stand for? 
Typical examples for NICE{n) are formulae such as 


1 


n -|- 1 

{Catalan numbers] cf. jl78l Ex. 6.19]) or 


2n 

n 


(3i + 1)! 


n—1 

n (n + i)\ 


(2.3) 


(2.4) 


^For more information on this determinant see Theorems |21 and ^ in this section and ITTO Sec¬ 
tions 2.2, 2.3 and 2.5]. 

^This determinant evaluation follows easily from the matrix factorisation 

(S(i + J,i))i<,,,<„ = ((-l)‘tV(f.'!(i - *)!)), <,K„ . (*').«,,<„ ■ 

application of ITTM Theorem 26, (3.14)] to the first determinant, and application of the Vandermonde 
determinant evaluation to the second. 
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(the number of n x n alternating sign matrices and several other combinatorial objects; 
cf. [2H1)- Let us be more precise. 


“Definition”. The symbol NICE{n) is a formula of the type 


C ■ Rat(n) 


n(ajn + bi)\ 
{ciU + dj)!’ 

2=i ^ ^ 


(2.5) 


where Rat(n) is a rational function in n, and where ai, Ci & T, for i = 1, 2,..., fc, Z 
denoting the set of integers. The parameters can be arbitrary real or complex 

numbers. {If necessary, {uiU + hi)\ has to be interpreted as T{ain + bi + 1), where r(a;) 
is the Euler gamma function, and similarly for {ciU + di)\.) 


Clearly, the formulae (Q and (12.411 £t this “Dehnition”.^ 

If one is working on a particular problem, how can one recognise that one is looking at 
a sequence of numbers given by NICE{n)l The key observation is that, if we factorise 
{an + b)! into its prime factors, where a and b are integers, then, as n runs through 
the positive integers, the numbers {an + 6)! explode quickly, whereas the prime factors 
occurring in the factorisation will grow only moderately, more precisely, they will grow 
roughly linearly. Thus, if we encounter a sequence the prime factorisation of which has 
this property, we can be sure that there is a formula NICE{n) for this sequence. Even 
better, as I explain in Appendix A of |l()9j . the program Rate^ will (normally®) be able 
to guess the formula. 

To illustrate this, let us look at a particular example. Let us suppose that the first 
few values of our sequence are the following: 


1,2,5,14,42,132,429,1430,4862,16796, 58786, 208012, 742900, 2674440, 

9694845, 35357670,129644790,477638700,1767263190,6564120420. 

The prime factorisation of the second-to-last number is (we are using Mathematica 
here) 

In[l]:= Factorinteger[477638700] 


^The writing NICE{n) is borrowed from Doron Zeilberger [inni Recitation III]. The technical term 
for a formula of the type (12.511 is “hypergeometric term”, see [mi Sec. 3.2], whereas, most often, the 
colloquial terms “closed form” or “nice formula” are used for it, see dnsi Recitation II]. More recently, 
some authors call sequences given by formulae of that type sequences of “round” numbers, see [TT71 
Sec. 6]. 

%ate is available from http://igd.univ-lyonl.fr/~kratt It is based on a rather simple algo¬ 
rithm which involves rational interpolation. In contrast to what I read, with great surprise, in gni, 
the explanations of how Rate works in Appendix A of unsi can be read and understood without any 
knowledge about determinants and, in particular, without any knowledge of the fifty or so pages that 
precede Appendix A in unni. 

®Rate will always be able to guess a formula of the type (1^ if there are enough initial terms of 
the sequence available. However, there is a larger class of sequences which have the property that 
the size of the primes in the prime factorisation of the terms of the sequence grows only slowly with 
n. These are sequences given by formulae containing “Abelian” factors, such as n". Unfortunately, 
Rate does not know how to handle such factors. Recently, Rubey HHBl proposed an algorithm for 
covering Abelian factors as well. His implementation Guess is written in Axiom and is available at 
http://www.mat.univie.ac.at/~rubey/martin.html 
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Out[l]= {{2, 2}, {3, 1}, {5, 2}, {7, 1}, {11, 1}, {23, 1}, {29, 1}, {31, l}} 

whereas the prime factorisations of the next-to-last and the last nnmber in this sequence 
are 

In[2]:= Factorinteger[1767263190] 

0ut[2]= {{2, 1}, {3, 1}, {5, 1}, {7, 1}, {11, 1}, {23, 1}, {29, 1}, {31, 1}, 

> { 37 , 1 }} 

In[3]:= Factorinteger[6564120420] 

□ut[3]= {{2, 2}, {3, 1}, {5, 1}, {11, 1}, {13, 1}, {23, 1}, {29, 1}, {31, 1}, 

> { 37 , 1 }} 

(To decipher this for the reader unfamiliar with Mathematica: the prime factorisation 
of the last number is 2^3^5^11^13^23^29^31^37^.) One observes, hrst of all, that the 
occurring prime factors are rather small in comparison to the numbers of which they 
are factors, and, second, that the size of the prime factors grows only very slowly (from 
31 to 37). Thus, we can be sure that there is a “nice” formula NICE{n) for this 
sequence. Indeed, Rate needs only the first five members of the sequence to come up 
with a guess for NICE{n): 

In [4] : = «rate .m 
In[5]:= Rate[1,2,5,14,42] 
iO 1 

4 Gamma[- + iO] 

2 

□ut[5]=- 

Sqrt[Pi] Gamma[2 + iO] 

As the reader will have guessed. Rate uses the parameter zq instead of n. In fact, the 
formula is a fancy way to write that is, we were looking at the sequence of 

Catalan numbers (EH- 

To see the sharp contrast, here are the first few terms of another sequence: 

1,2,9,272,589185. 

(Also these are combinatorial numbers. They count the perfect matchings of the n- 
dimensional hypercube; cf. |14fil Problem 19].) Let us factorise the last two numbers: 

In[6]:= Factorinteger[272] 

□ut[6]= {{2, 4}, {17, 1}} 

In[7]:= Factorinteger [589185] 

0ut[7]= {{3, 2}, {5, 1}, {13093, l}} 

The presence of the big prime factor 13093 in the last factorisation is a sure sign that 
we cannot expect a formula NICE{n) as described in the “Definition” for this sequence 
of numbers. (There may well be a simple formula of a different kind. It is not very 
likely, though. In any case, such a formula has not been found up to this date.) 
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Now, that I have sufficiently explained all the ingredients in the “prototype theorem” 
Theorem ^ I can explain why theorems of this form are so attractive (at least to me): 
the objects (i.e., the permutations, partitions, paths, trees, etc.) that it deals with are 
usually very simple to explain, the statement is very simple and can be understood by 
anybody, the result NICE{n) has a very elegant form, and yet, very often it is not 
easy at all to give a proof (not to mention a true explanation why such an elegant result 
occurs.) 

Here are two examples. They concern rhombus tilings, by which I mean tilings of a 
region by rhombi with side lengths 1 and angles of 60° and 120°. The hrst one is a one 
century old theorem due to MacMahon |ldl[ Sec. 429, q —>■ 1] proof in Sec. 494].® 


Theorem 2. The number of rhombus tilings of a hexagon with side lengths a, b, c, a, b, c 
whose angles are 120° (see Figure\^a for an example of such a hexagon, and Figure\^b 
for an example of a rhombus tiling) is equal to 


n 


{a + b + i — 1) \ {i — 1)! 
{a + i — 1) \ {b + i — 1)! 


( 2 . 6 ) 

□ 



a. A hexagon with sides a, b, c, a, b, c, 
where a = 3, 6 = 4, c = 5 


b. A rhombus tiling of a hexagon 
with sides a, b, c, a, 6, c 


Figure 1. 


®To be correct, MacMahon did not know anything about rhombus tilings, they did not exist in 
enumerative combinatorics at the time. The objects that he considered were plane partitions. However, 
there is a very simple bijection between plane partitions contained in an a x 6 x c box and rhombus 
tilings of a hexagon with side lengths a, b, c, a, b, c, as explained for example in EH 
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The second one is more recent, and is due to Ciucu, Eisenkolbl, Zare and the author 
m Theorem 1], 

Theorem 3. If a, b, c have the same parity, then the number of lozenge tilings of a 
hexagon with side lengths a,b + m,c,a + m,b,c + m, with an equilateral triangle of side 
length m removed from its centre {see Figure\^for an example) is given by 


H(a 

X 


m 


H(6 + m)H(c + m)H(a + & + c + m) H(m + \^±^])R{m+ 
6 + m) H(a + c + m) H(6 + c + m) H(^ + m)H(^^ + ; 


Rja + b + c + m) H(m+ [g±|±£])H(m+ [^J) 
i + c + m) H(6 + c + m) H(^ + m) + m) H(^ + rr 

H(f+ [ID H(f + [ID H(f + [fD H(f + [|J) H(f + [|J) H( 

7n\2 g+fe+'ffl- ^2 jj{ a+c+m \2 jj{ fc+c+m \2 


m) 


X 


H(f) 


TT/ m 




- 2 

I” y~f 


g+c+m ^2 fc+c+m ^2 


+ [fj) 


where 


H(n) 


re;ir('=+i) 
nrJ nk + 1 ) 


for n an integer, 
for n a half-integer. 


( 2 . 8 ) 

□ 


(There is a similar theorem if the parities of a, b, c should not be the same, see |4()1 
Theorem 2]. Together, the two theorems generalise MacMahon’s Theorem |21^) 

The reader should notice that the right-hand side of (ESD is indeed of the form 
NICE{a), while the right-hand side of ()2.7I1 is of the form NICE{m/2). 

Where is the connexion to determinants? As it turns out, these two theorems are in 
fact determinant evaluation theorems. More precisely. Theorem |21 is equivalent to the 
following theorem. 

'^Bijective proofs of Theorem |3 which “explain” the “nice” formula are known [lOfiLll()ll| . I do not 
ask for a bijective proof of Theorem |21 because I consider the task of finding one as daunting. 
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Theorem 4. 

det f f ^ ^ ~ ~ 1)! 

i<*j<c \ \a — i + j) ) (a + i — 1)! (6 + i — 1)! 


(2.9) 

□ 


(The reader should notice that this is exactly ()2.1|1 with n replaced by c.) On the 
other hand, Theorem IHl is equivalent to the theorem below.® 


Theorem 5. If m is even, the determinant 

/ fb + c + m 


det 

l<ij<a+m 


b — i + j 

b-\-c 
2 


1 < i < a 


\ 


a+l<i<a+m 


is equal to (EH). 


( 2 . 10 ) 

□ 


The link between rhombus tilings (and equivalent objects such as plane partitions, 
semistandard tableaux, etc.) and determinants which explains the above two equiva¬ 
lence statements is non-intersecting lattice paths.^ The latter are families of paths in 
a lattice with the property that no two paths in the family have a point in common. 
Indeed, rhombus tilings are (usually) in bijection with families of non-intersecting paths 
in the integer lattice which consist of unit horizontal and vertical steps. (Figure El 
illustrates the bijection for the rhombus tilings which appear in Theorem El in an exam¬ 
ple. In that bijection, all horizontal steps of the paths are in the positive direction, and 
all vertical steps are in the negative direction. See the explanations that accompany uni 
Figure 8] for a detailed description. Since, as I explained, TheoremElessentially is a gen¬ 
eralisation of Theorem El this gives also an idea for the bijection for the rhombus tilings 
which appear in the latter theorem. For other instances of bijections between rhombus 
tilings and non-intersecting lattice paths see [¥n 1^iKKi iKHi IK71 111211141] . 1. In 
the case that the starting points and the end points of the lattice paths are fixed, the 
following many-author-theorem applies. 

®To be correct, this is a little bit oversimplified. The truth is that equivalence holds only if m is 
even. An additional argument is necessary for proving the result for the case that m is odd. We refer 
the reader who is interested in these details to @ni Sec. 2]. 

®There exists in fact a second link between rhombus tilings and determinants which is not less 
interesting or less important. It is a well-known fact that rhombus tilings are in bijection with perfect 
matchings of certain hexagonal graphs. (See for example fTTHl Figures 13 and 14].) In view of this 
fact, this second link is given by Kasteleyn’s theorem m saying that the number of perfect matchings 
of a planar graph is given by the Pfafhan of a slight perturbation of the adjacency matrix of the graph. 
See | llfi| for an exposition of Kasteleyn’s result, including historical notes, and for adaptations taking 
symmetries of the graph into account. 

^^This result was discovered and rediscovered several times. In a probabilistic form, it occurs for the 
first time in work by Karlin and McGregor In matroid theory, it is discovered in its discrete 

form by Lindstrom Lemma 1]. Then, in the 1980s the theorem is rediscovered at about the same 
time in three different communities, not knowing from each other at the time: in statistical physics by 
Fisher Sec. 5.3] in order to apply it to the analysis of vicious walkers as a model of wetting and 
melting, in combinatorial chemistry by John and Sachs inn and Gronau, Just, Schade, Scheffler and 
Wojciechowski m in order to compute Pauling’s bond order in benzenoid hydrocarbon molecules, 
and in enumerative combinatorics by Gessel and Viennot cm El in order to count tableaux and plane 
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a. A lozenge tiling of the cored hexagon in Figure I^J b. The corresponding path family 



Figure 3. 

Theorem 6 (Karlin-McGregor, Lindstrom, Gessel-Viennot, Fisher, 
John-Saghs, Gronau-Just-Sghade-Sgheeeler-Wojgieghowski). Let Ai , A2 , 
^ An and Ei, E2, ..., En be lattice points such that for i < j and k < I any lattice path 
between Ai and Ei has a common point with any lattice path between Aj and Ej. . Then 
the number of all families (Pi, P 2 ,..., Pn) of non-intersecting lattice paths, Pi running 
from Ai to Ei , i = 1 , 2 ,. .. ,n, is given by 

where P{A —>■ E) denotes the number of all lattice paths from A to E. □ 

It goes beyond the scope of this article to include the proof of this theorem here. 
However, I cannot help telling that it is an extremely beautiful and simple proof that 

partitions. Since only Lindstrom, and then Gessel and Viennot state the result in its most general 
form (not reproduced here), I call this theorem most often the “Lindstrom-Gessel-Viennot theorem.” 
It must be also mentioned that the so-called “Slater determinant” in quantum mechanics (cf. |170| and 
irzn Ch. 11]) may qualify as an “ancestor” of the Lindstrom-Gessel-Viennot determinant. 




































ADVANCED DETERMINANT CALCULUS: A COMPLEMENT 


9 


every mathematician should have seen once, even if (s)he does not have any use for it 
in her/his own research. I refer the reader to |73l 17^ ri79j . 

Now the origin of the determinants becomes evident. In particular, since, for rhom¬ 
bus tilings, we have to deal with lattice paths in the integer lattice consisting of unit 
horizontal and vertical steps, and since the number of such lattice paths which connect 
two lattice points is given by a binomial coefficient, we see that the enumeration of 
rhombus tilings must be a rich source for binomial determinants. This is indeed the 
case, and there are several instances in which such determinants can be evaluated in 
the form (see j!?KlHHlITTl H!?!liH IHHlIHTHIH71 Fin ITTl |l()9j and Sectional). Often 

the evaluation part is highly non-trivial. 

The evaluation of the determinant (im is not very difficult (see [TTM Sections 2.2, 
2.3, 2.5] for 3 different ways to evaluate it). On the other hand, the evaluation of the 
determinant (ITTUl) requires some effort (see m Sec. 7]). 

To conclude this section, I state another determinant evaluation, to which I shall 
come back later. Its origin lies as well in the enumeration of rhombus tilings and plane 
partitions (see |1()51 Theorem 10] and jl2l Theorem 2.1]). 

Theorem 7. For any complex numbers x and y there holds 

det ( (x + v + , + j-l)\ \ 

o<i,j<n-i \ (x + 2i — j)\ [y -\- 2j — i) \ J 

_'^j^i\{x + y + i-l)\{2x + y + 2i)i {x + 2y + 2i)i 

{x + 2i)\{y + 2i)\ ’ I • 1 

where the shifted factorials or Pochhammer symbols (a)^ are defined by {a)k '■= a{a -|- 
1) ■ ■ ■ (a -|- fc — 1), k > 1, and (a)o := 1. {In this formula, a factorial ml has to be 
interpreted as r(m -|- 1) if m is not a non-negative integer.) 


3. A determinant from number theory 


However, determinants do not only arise in combinatorics, they also arise in other 
helds. In this section, I want to present a determinant which arose in number theory, 
explain in some detail its origin, and then outline the steps which led to its evaluation, 
thereby giving the reader an opportunity to look “behind the scenes” while one tries to 
make the determinant evaluation methods described in uns] work. 

The story begins with the following two series expansions for tt. The hrst one is due 
to Bill Gosper UTI, 


TT = 


oo 


E 

n=0 


50n-6 


(3.1) 


and was used by Fabrice Bellard |2()1 hie pil.c] to hnd an algorithm for computing 
the n-th decimal of vr without computing the earlier ones, thus improving an earlier 
algorithm due to Simon Plouffe |145j . The second one. 


1 

740025 


3P(n) 

\n=l \2n)^ 


20379280j 


(3.2) 
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where 

P{n) = -885673181n^ + 3125347237n^ - 2942969225n^ 

+ 1031962795n2 - 196882274n + 10996648, 

is due to Fabrice Bellard [201, and was used by him in his world record setting compu¬ 
tation of the 1000 billionth binary digit of vr, being based on the algorithm in |19j . 

Going beyond that, my co-authors from j3], Gert Almkvist and Joakim Petersson, 
asked themselves the following question: 

Are there more expansions of the type 


TT = 




S(n) 


n=0 \pn/ 

where S{n) is some polynomial in n {depending on m,p, a) ? 

How can one go about to get some intuition about this question? One chooses some 
specihc m,p,a, goes to the computer, computes 




n 


/ j lmn\ „ 

— n { O 

n=0 \pn/ 

to many, many digits for fc = 0,1, 2,..., puts 

7r,p(0),p(l),p(2),... 

into the LLL-algorithm (which comes, for example, with the Maple computer algebra 
package), and one waits whether the algorithm comes up with an integral linear com¬ 
bination of 7r,p(0),p(l),p(2),... Indeed, Table [T] shows the parameter values, where 
the LLL-algorithm gave a result. 

For example, it found 

Iv^ S{n) 


TT = 


n=0 


where 


r = 3®5^7^11^13^ 


and 


S{n) = -869897157255 - 3524219363487888n + 112466777263118189n2 
- 1242789726208374386/1^ + 6693196178751930680n^ - 19768094496651298112n® 
+ 32808347163463348736/1® - 28892659596072587264n^ + 10530503748472012800n®, 


llpor readers unfamiliar with the LLL-algorithm: in this particular application, it takes as an 
input rational numbers ri,r 2 ,...,rm (which, in our case, will be the numbers 1 and the rational 
approximations of tt, p(0), p(l), . .. which we computed), and, if successful, outputs small integers 
Cl, C 2 , ..., Cm such that ciri -I-C 2 r 2 -|- • • • -t-Cmim is very small. Thus, if was a good approximation for 
the real number i = 1, 2,..., m, one can expect that actually cia:i -|-C 2 a ;2 -I-- • • + Cma:m = 0. See |12dl 
Sec. 1, in particular the last paragraph] and usi Ch. 2] for the description of and more information 
on this important algorithm. In particular, also here, the output of the algorithm (if there is) is just 
a (very guided) guess. Thus, a proof is still needed, although the probability that the guess is wrong 
is infinitesimal. As a matter of fact, it is very likely that Bellard had no proof of his formula ... 






ADVANCED DETERMINANT CALCULUS: A COMPLEMENT 


11 


m 

P 

a 

deg(E) 

3 

1 

2 

1 

7 

2 

2 

5 

8 

4 

-4 

4 

10 

4 

4 

8 

12 

4 

-4 

8 

16 

8 

16 

8 

24 

12 

-64 

12 

32 

16 

256 

16 

40 

20 

_45 

20 

48 

24 

46 

24 

56 

28 

-47 

28 

64 

32 

48 

32 

72 

36 

-49 

36 

80 

40 

410 

40 


(Gosper) 

(Bellard) 


Table 1. 


and 


where 


and 


TT 


1 

r 


OO 


E 

n =0 


S{n) 

(? 3256 ^ 


r = 2^3^^5^7hlh3^17^19^23‘^29^3l‘^ 


S{n) = - 2062111884756347479085709280875 

+ 1505491740302839023753569717261882091900n 

- 112401149404087658213839386716211975291975n2 
+ 3257881651942682891818557726225840674110002n^ 
-51677309510890630500607898599463036267961280n^ 

+ 517337977987354819322786909541179043148522720n^ 

- 3526396494329560718758086392841258152390245120n® 

+ 171145766235995166227501216110074805943799363584n^ 

- 60739416613228219940886539658145904402068029440n® 

+ 159935882563435860391195903248596461569183580160n^ 
-313951952615028230229958218839819183812205608960n^° 
+ 457341091673257198565533286493831205566468325376n^^ 
-486846784774707448105420279985074159657397780480n^2 
+367314505118245777241612044490633887668208926720n^^ 

- 185647326591648164598342857319777582801297080320n^^ 
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+ 56224688035707015687999128994324690418467340288n^^ 
- 7687255778816557786073977795149360408612044800n^® 


Of course, there could be many more. 

If one looks more closely at Table d then, if one disregards the hrst, second and 
fourth line, one cannot escape to notice a pattern: apparently, for each k = 1 , 2 ,..., 
there is a formula 

_ Sk{n) 

^ ~ 2-^ (^kn\/_A\kn' 
n=0 \4kn/ ^ 

where Sk{n) is some polynomial in n of degree Ak. 

In order to make progress on this observation, we have to hrst see how one can prove 
such an identity, once it is found. In fact, this is not difficult at all. To illustrate the 
idea, let us go through a proof of Gosper’s identity (EH). 

The beta integral evaluation (cf. [T^ Theorem 1.1.4]) gives 

1 

y^ = (3u + l) / x^’^{l-xYdx. 

(n) 

Hence the right hand side of the formula will be 


„1 oo 

/ -6)(3n + l) 

n=0 


x‘^{l — x) 


dx. 


We have 


— 6)(3n + Vjy"^ 


2(56i/2 + 97?/ - 3) 


n=0 


(1-2/)^ 

Thus, if substituted, we obtain 

2Sx^ - 56x5 ^ 28x^ - 97x3 ^ 97^2 _ g 


(3.3) 


RHS = 


(x3 _ ^2 + 2)3 


dx 


4 x(x — l)(x 3 — 28 x^ + 9 x + 8) ^ ^ 

^^1 + 4 arctan X - 1 
(X'^ — x2 + 2)2 


= TT. 


(3.4) 


(Clearly, both (I3.3|l and dH are routine calculations, and therefore we did not do it 
by hand, but let them be worked out by Maple.) 

Now let us hx fc > 1. We apply the same procedure to '222=0 
where Sk{n) is (hopefully) some (unknown) polynomial in n. The beta integral evalu¬ 
ation gives 

1 


/8kn\ 

\4fcn/ 


= {8kn + l) / x^'”^(l-x)^'”"dx. 


Hence, if Sk{n) should have degree d in n. 

Skin) 




(_A)kn L A_ 
=0 \iknj\ n=< 


'^^{Skn -|- l)S'j 


k[n) 


Pkix) 


lo (x‘^^(l — x)^^ — (—4)*^)'^"''^ 


x4fe(i _ 

~TW 

dx, 


4k 


dx 


(3.5) 
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where Pk{x) is some polynomial in x. For convenience, let us write P as a short-hand 
for Pfc. Let Q{x) := x^^{l — xY^ — (—4)^. Now we make the wild assumption that 

f Pix) R(x) 

/ dx = . -|- 2arctanfa;) -|- 2arctanfa; — 1), 

J Q[xY+^ Q[xY+^ 

for some polynomial R{x) with P(0) = P(l) = 0. Then the original sum would indeed 
be equal to tt. The last equality is equivalent to 


P 


R' 


Qd+2 Qd+l 




+ 


or 


QR! - (d + 1)Q'P = P - 2Q'^+2 
In our examples, we observed that 


-I- 1 — 2x -I- 2 

1 1 
+ 


xi -1- 1 — 2x + 2 


R{x) = (2x — l)P(x(l — x)) 
for a polynomial R. So, let us make the substitution 

t = x(l — x). 


Then, after some simplihcation, the above differential equation becomes 

ap p)4A:+2 

-(1 - 4t)g^ + (2g + Ak{Ak + 1)(1 - - P + 2(3 - 2t )= 0, (3.6) 

where Q{t) = - (-4)*=. 

Now, writing N{k) = 4k{4k -l- 1), we make the Ansatz 

N{k)-l 

Rit) = 

3=^ 

Ak 

Sk{n) = J2<N{k)+j)n^. 

j=0 

(The reader should recall that Sk{n) dehnes Pkit) = P{t) through (I3.5j) .l Comparing 
coefficients of powers of t on both sides of (EnD, we get a system of N{k) -|- 4k linear 
equations for the unknowns a(l), a(2),..., a{N{k) -|- 4k). 

Hence: If the determinant of this system of linear equations is non-zero, then there 
does indeed exist a representation 

Rk{n) 

n=0 \4fcn/ ^ ' 


To see whether we could indeed hope for the determinant to be non-zero, we went 
again to the computer and looked at the values of the determinant in some small 
instances. (Obviously, we do not want to do this by hand, since for k = 1 the matrix is 
already a 24x 24 matrix!) So, let us program the matrix. (We shall see the mathematical 
dehnition of the matrix in just a moment, see (EHl)-^^) 

tell the truth, this is the form of the matrix after some simplifications have already been carried 
out. (In particular, we are looking at a matrix which is slightly smaller than the original one.) See 13 
beginning of Section 4] for these details. There, the matrix in 13.811 is called M'". 
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In[8]:= a[k_,j_]:=Module[{Var=j/(4k)}, 

(-l)^(Var-l)*8k(4k+l)(-4)^(k*(Var+l))* 

Product[4k*l-l,{l,1,4k-Var}]^Product[4k*l+l,{l,1,Var-l}] 

] 


In[9]:= A[k_, i_, j_] := Module[{Var}, 

Var = {Floor[(i - 2)/(4*k - 1)], 

Floor[(j - l)/(4*k)], Mod[i - 2, 4*k - 1], 
Mod[j - 1, 4*k]}; 

If[i == 1, 


If [Mod[j, 4*k] === 0, a[k, j] , 0], 

If[Var[[l]] - Var[[2]] == 0, 

Switch[Var[[3]] - Var[[4]], 0, fl[k, Var[[3]]+1, j], -1, 
f0[k, Var[[3]]+1, j] , 0], 

If[Var[[1]] - Var[[2]] == 1, 

Switch[Var[[3]] - Var[[4]], 0, gl[k, Var[[3]]+1, j], -1, 
g0[k, Var[[3]]+1, j] , 0], 0]]]] 

In[10]:= A[k_] :=Table[A[k, i, j], {i, 1, 16*k^2}, {j, 1, 16*k^2}] 
In[ll]:= f0[k_, t_, j.] := j*(-4)^k; 

fl[k_, t_, j.] := -(2 + 4*j)*(-4)^k; 

g0[k_, t_, j.] := (4=i'k=i'(4=i'k + 1) - j); 

gl [k_, t_, j_] := (-4*4*k*(4*k + 1) + 2 + 4*j) 


We shall not try to digest this at this point. Let us accept the program as a black 
box, and let us compute the determinant for k = 2. 

In [12] := Det[A[2]] 

Out[12]= -601576375580370166777074138698518196031142518971568946712\ 

> 2204136674781038302774231725971306459064075121023092662279814\ 

> 015195545600000000000 

Magnihcent! This is certainly not zero. However, what are we going to do with this 
gigantic number? Remembering our discussion about “nice” numbers and “nice” for¬ 
mulae in the preceding section, let us factorise it in its prime factors. 

In[13]:= Factorinteger[%] 

0ut[13]= {{-1, 1}, {2, 325}, {3, 39}, {5, 11}, {7, 11}, {11, 3}, {13, 2}} 

1 would say that this is sensational: a number with 139 digits, and the biggest prime 
factor is 13! As a matter of fact, this is not just a rare exception. Table El shows the 
factorisations of the hrst hve determinants. (We could not go further because of the 
exploding size of the matrix of which the determinant is taken.) 
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k 

1 

2 

3 

4 

5 


det(y4(A;)) 

259355671 

_2325 3395ll7llii3i32 

27723I465287I777I773I877479323I 

-2I9133III558738772I7322772479723529231I 

22932320253067697729732777287929239296375372 


Table 2. 


Thus, these experimental results make us sure that there must be a “nice” formula 
for the determinant. Indeed, we prove in that^^ 

det{A{k)) = + 1)!^^ JJ (3.7) 

j=i 

Hence the desired theorem follows. 

Theorem 8. For all k > 1 there is a formula 

Sk{n) 

n=0 \4kn/ ^ ^ 

where Sk(n) is a polynomial in n of degree Ak with rational coefficients. The polynomial 
Sk(n) can be found by solving the previously described system of linear eguations. 

I must admit that we were extremely lucky that it was indeed possible to evaluate the 
determinant explicitly. To recall, “all” we needed to prove our theorem (Theorem IHj) 
was to show that the determinant was non-zero. To be honest, I would not have the 
slightest idea how to do this here without hnding the exact value of the determinant. 

Now, after all this somewhat “dry” discussion, let me present the determinant. We 


had to determine the determinant of the 

X matrix 



/O ... 0 * 

0...0* 0...0* ... 



* 

O 

o 



El 

0 

0 ... 



0 



Gi 

F2 

0 ... 



0 



0 

G2 

Es 






0 

0 

Gs 




. (3.8) 






^4fc-l 

0 





0 

0 

G/ik-i 

F^k 



0 



0 

0 

Gik / 


^^Strictly speaking, this is not a 

formula NICE(k) according to my 

“Definition’ 

in the preceding 


section, because of the presence of the “Abelian” factors and (4k + 1)!^^, see Footnote El 

Nevertheless, the reader will certainly admit that this is a nice and closed formula. 
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where the £-th non-zero entry in the hrst row (these are marked by *) is 






n(4«-i) 


2 = 1 


n-1 


]^(4ifc -h 1) 


y2=l 




and where each block Ft and Gt is a (4fc —1) x {4k) matrix (that is, these are rectangular 
blocks!) with non-zero entries only on its (two) main diagonals, 

//l(4(t-l)fc+l) /o(4(t-l)fc+2) 0 ... \ 

0 /i(4(t-l)fc+2) /o(4(4-l)fc+3) 0 


Ft = 


\ 


0 /i(41A:-2) /o(41A:-1) 0 

0 /i(4tfc-l) fo(4:tk) ) 


and 


/Sl(4(t-l)fc+l) go(4(t-l)fc+2) 0 

0 gi(4(t-l)fc+2) go(4(t-l)fc+3) 0 


Gt = 


\ 


0 gi{Atk—2) 5'o(4tfc—1) 0 

\ 0 gi{uk-l) go(itk) ) 

We have almost worked our way through the dehnition of the determinant. The only 
missing piece is the dehnition of the functions /o, /i, ( 70 , 9\ in the blocks Ft and Gt- Here 
it is: 


/o(J)=J(-4)^ 

/l(J) = -(4J + 2)(-4)^ 

77o(j) = (At(fc)-i), 

^i(j) = -(4iV(fc)-4j-2), (3.9) 

where, as before, we write At(/c) = 4k{4k + 1) for short. 


4. The evaluation of the determinant 

I now describe how the determinant of (HHI) was evaluated by applying to it the 
methods described in [inn] . To make this section as self-contained as possible, for each 
of them I briehy recall how it works before putting it into action. 

“Method” 0: Do row and column operations until the determinant reduces to some¬ 
thing manageable. 

In fact, at a hrst glance, this does not look too bad. Our matrix (EH), of which we 
want to compute the determinant and show that it is non-zero, is a very sparse matrix. 
Moreover, it looks almost like a two-diagonal matrix. It seems that one should be able 
to do a few row and column manipulations and thus reduce the matrix to a matrix of 
a simpler form of which we can evaluate the determinant. 

Well, we tried that. Unfortunately, the above impression is deceiving. First of all, 
the diagonals of the blocks do not really ht together to form diagonals which run from 
one end of the matrix to the other. Second, there remains still the hrst row which does 
not ht the pattern of the rest of the matrix. So, whatever we did, we ended up nowhere. 
Maybe we should try something more sophisticated ... 
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Method 1 jl(J9L Sec. 2.6]: LU-factorisation. Suppose we are given a family of matrices 
y4(l), ^4(2), y4(3),... of which we want to compute the determinants. Suppose further 
that we can write 

A{k) ■ U{k) = L{k), 

where U{k) is an upper triangular matrix with Is on the diagonal, and where L{k) is a 
lower triangular matrix. Then, clearly, 

det(A(/c)) = product of the diagonal entries of L{k). 

But how do we find U{k) and L{k)l We go to the computer, crank out U{k) and 
L{k) for k = 1, 2, 3,..., until we are able to make a guess. Afterwards we prove the 
guess by proving the corresponding identities. 

Well, we programmed that, we stared at the output on the computer screen, but we 
could not make any sense of it. 

Method 2 |1()91 Sec. 2.3]: Condensation. This is based on a determinant formula due 
to Jacobi (see |2Hl Ch. 4] and |1()()1 Sec. 3]). Let A be an n x n matrix. Let 
denote the submatrix of A in which rows ii,i 2 , ■ ■ ■ Ci and columns ji, j 2 ,..., ji are 
omitted. Then there holds 

det A ■ det A}’” = det Aj ■ det A" — det A" • det A^. (4.1) 

If we consider a family of matrices A(l), A(2),..., and if all the consecutive minors of 
A(n) belong to the same family, then this allows one to give an inductive proof of a 
conjectured determinant evaluation for A(n). 

Let me illustrate this by reproducing Amdeberhan’s condensation proof [S] of ()2.11|1 . 
Let Mn{x,y) denote the determinant in ()2.11|1 . Then we have 

{Mn{x,y))l = Mn-i{x,y), 

{Mn{x,y))l = Mn-i{x + l,y + 1), 

{Mn{x, y))l = Mn-i{x -l,y+ 2), 

{Mn{x, y))'l = Mn-i{x + 2,y-l), 

{Mr,{x,y)Yf^ = Mn-2{x + l,y + l). (4.2) 

Thus, we know that Equation (ED is satisfied with A replaced by Mn{x,y), where the 
minors appearing in (ED are given by (03). This can be interpreted as a recurrence for 
the sequence [Mn{x,y))^.^^. Indeed, given Mo{x,y) and Mi{x,y), the equation (14.Ij) 
determines Mn{x,y) uniquely for all n > 0 (given that Mn{x,y) never vanishes). Thus, 
since the right-hand side of (imi) is indeed never zero, for the proof of it suffices 

to check (ITTTl) for n = 0 and n = 1, and that the right-hand side of (imu also satisfies 
(ED, all of which is a routine task. 

Now, a short glance at the definition of our matrix (EHD will convince us quickly that 
application of this method to it is rather hopeless. For example, omission of the first 
row already brings us outside of our family of matrices. So, also this method is not 
much help to solve our problem, which is really a pity, because it is the most painless 

of all ... 
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Method 3 |1(J91 Sec. 2.4]: Identification of factors. In order to sketch the idea, let us 
quickly go through a (standard) proof of the Vandermonde determinant evaluation, 

n (4.3) 

l< 2 <j<n 

Proof. If = Xi^ with ii ^ i 2 , then the Vandermonde determinant (I4.d|l certainly 
vanishes because in that case two rows of the determinant are identical. Hence, — 
Xjj) divides the determinant as a polynomial in the Xj’s. But that means that the 
complete product ni<i<j<n(^i “ (which is exactly the right-hand side of ()4.dj) i 
must divide the determinant. 

On the other hand, the determinant is a polynomial in the Xj’s of degree at most 
( 2 ). Combined with the previous observation, this implies that the determinant equals 
the right-hand side product times, possibly, some constant. To compute the constant, 
compare coefficients of X^X\ ■ ■ ■ Xtf~^ on both sides of (031). This completes the proof 
of (j4.3|) . □ 

At this point, let us extract the essence of this proof. The basic steps are: 

(51) Identification of factors 

(52) Determination of degree bound 

(53) Computation of the multiplicative constant. 

As I report in uni, this turns out to be an extremely powerful method which has 
numerous applications. To given an idea of the flavour of the method, I show a few 
steps when it is applied to the determinant in ()2.11|1 (ignoring the fact that we have 
already found a very simple proof of its evaluation; see una proof of Theorem 10] for 
the complete proof using the “identihcation of factors” method). 

To get started, we have to transform the assertion (ITTTD into an assertion about 
polynomials. This is easily done, we just have to factor 

{x + y + i — V)\/{x + ‘2i)\/{y + 2n — i — 2)\ 

out of the f-th row of the determinant. If we subsequently cancel common factors on 
both sides of dTTTl) . we arrive at the equivalent assertion 

det ((x + 2 / + i)j {x + 2i- j + 1)^ {y + 2j-i + l) 2 n- 2 j- 2 ) 

0 < 2 ,j<n—1 ' 

n—1 

= (*! ( 2 / + -|- l)n-i-i {2.x + y + 2i)i {x + 2y + 2i)fj , (4.4) 

i=0 

where, as before, {a)k is the standard notation for shifted factorials (Pochhammer 
symbols) explained in the statement of Theorem [Tj 

In order to apply the same idea as in the above evaluation of the Vandermonde 
determinant, as a hrst step we have to show that the right-hand side of (031) divides 
the determinant on the left-hand side as a polynomial in x and y. For example, we 
would have to prove that {x+2y+2n—2) (actually, (x-f2j/-|-2n—2)L(”+^F3J, we will come 
to that in a moment) divides the determinant. Equivalently, if we set x = —2y — 2n + 2 
in the determinant, then it should vanish. How could we prove that? Well, if it vanishes 
then there must be a linear combination of the columns, or of the rows, that vanishes. 
Equivalently, for x = —2y — 2n + 2 we hnd a vector in the kernel of the matrix in 
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(lOI), respectively of its transpose. More generally (and this addresses the fact that we 
actually want to prove that (x + 2?/ + 2n — 2)L("-+^h3J divides the determinant): 

For proving that {x + 2y + 2n — 2)^ divides the determinant, we find E linear 
independent veetors in the kernel. 

(For a formal justihcation that this does indeed suffice, see Section 2 of |l()7j . and in 
particular the Lemma in that section.) 

Okay, how is this done in practice? You go to your computer, crank out these vectors 
in the kernel, for n = 1, 2, 3,..., and try to make a guess what they are in general. To 
see how this works, let us do it in our example. First of all, we program the kernel of 
the matrix in (031) with X = —2y — 2n + 2 (again, we are using Mathematica here).^^ 

In[14]:= p=Pochhammer; 

ni[i_ , j-j] *P[y+S’i'j + l-i, 2*n-2*j-2] *p[x-j + l+2i, j] ; 

V[n_]:=(x=-2y-2n+2; 

Var=Sum[c[j]*Table [m [i,j,n],{i,0,n-l}],{j,0,n-l}]; 

Var=Solve[Var==Table [0,{n}] ,Table[c[i],{i,0,n-l}]] ; 

Factor[Table[c[i],{i,0,n-l}]/.Var]) 


What the computer gives is the following: 

In [15] := V[2] 

Out [15]= {{-2 c[l] , c[l]}} 

In [16] := V[3] 

Out [16]= {{-2 

In [17] := V[4] 

c[2] , 

-c [2] , c 

[2]}} 

Out[17]= {{-2 

In [18] := V[5] 

c[3] , 

-3 c[3] , 

0, c[3]}} 

Out[18]= {{-2 

In [19] := V[6] 

c[4] , 

-5 c [4] , 

-2 c[3] - c[4] , c[3] , c[4]}} 

Out [19]= {{-2 

In [20] := V[7] 

c[5] , 

-7 c[5] , 

-2 (c[4] + 2 c[5]), -c[4], c[4], c[5]}} 

Out [20]= {{-2 

c [6] , 

-9 c[6] , 

-2 c[5] - 9 c[6], -3 c[5] - c [6] , 0, 


> c[5], c[6]}} 

At this point, the computations become somewhat slow. So we should help our com¬ 
puter. Indeed, on the basis of what we have obtained so far, it is “obvious” that, 
somewhat unexpectedly, y does not appear in the result. Therefore we simply set y 
equal to some random number, and then the computer can go much further without 
any effort. 

In[21]:= y=101 


the program, V[n] represents the kernel, which is clearly a vector space. In the computer 
output, it is given in parametric form, the parameters being the c [i] ’s. 
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In [22] := V [8] 

Out [22]= {{-2 c[7], -11 c[7], -2 (c[6] + 8 c [7]) , -5 (c[6] + c [7]) , 

> -2 c[5] - c[6], c[5], c[6], c[7]}} 

In [23] := V[9] 

Out [23]= {{-2 c[8], -13 c[8], -2 c [7] - 25 c [8] , -7 (c[7] + 2 c [8]) , 

> -2 c[6] - 4 c[7] - c[8], -c[6], c [6] , c [7] , c[8]}} 

In [24] := V[10] 

Out [24]= {{-2 c[9], -15 c[9], -2 (c[8] + 18 c[9]), -3 (3 c [8] + 10 c [9]), 

> -2 c[7] - 9 c[8] - 6 c[9], -3 c [7] - c [8] , 0, c [7] , c [8] , c[9]}} 

In[25] := V[ll] 

Out [25]= {{-2 c[10], -17 c[10], -2 c [9] - 49 c[10], -11 (c[9] + 5 c[10]), 

> -2 (c[8] + 8 c[9] +10 c[10]), -5 c [8] - 5 c [9] - c[10], 

> -2 c[7] - c[8], c[7], c[8], c [9] , c[10]}} 

Let us extract some information out of these data. For convenience, we write 
for the matrix in (1131) in the sequel. For example, by just looking at the coefficients of 
c[n — 1] appearing in V[n], we extract that 

the vector (— 2 , 1 ) is in the kernel of M 2 , 

the vector (— 2 , — 1 , 1 ) is in the kernel of M 3 , 

the vector (—2, —3, 0,1) is in the kernel of M 4 , 

the vector (—2, —5, —1, 0,1) is in the kernel of M 5 , 

the vector (— 2 , —7, —4, 0, 0,1) is in the kernel of Mg, 

the vector (—2, —9, —9, —1, 0, 0,1) is in the kernel of My, 

the vector (—2, —11, —16, —5, 0, 0, 0,1) is in the kernel of Mg, 

the vector (—2, —13, —25, —14, —1, 0, 0, 0,1) is in the kernel of Mg, 

the vector (—2, —15, —36, —30, — 6 , 0, 0, 0, 0,1) is in the kernel of Mig, 

the vector (—2, —17, —49, —55, —20, —1, 0, 0, 0, 0,1) is in the kernel of Mu. 

Okay, now we have to make sense out of this. Our vectors in the kernel have the 
following structure: hrst, there are some negative numbers, then follow a few zeroes, 
and hnally there is a trailing 1. I believe that we do not have any problem to guess 
what the zeroeth^® or the hrst coordinate of our vector is. Since the second coordinates 
are always negatives of squares, there is also no problem there. What about the third 
coordinates? Starting with the vector for My, these are —1, —5, —14, —30, —55,.... I 
guess, rather than thinking hard, we should consult Rate (see Footnote il): 


^^The indexing convention in the matrix in of which the determinant is taken is that rows and 
columns are indexed by 0,l,...,n — 1. We keep this convention here. 
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In[26]:= Rate[-1,-5,-14,-30,-55] 

-(iO (1 + iO) (1+2 iO)) 

Out [26]= {-} 

6 

After replacing iO by n — 6 (as we should), this becomes —{n — 6 )(n — 5)(2n — ll)/ 6 . 
An interesting feature of this formula is that it also works well for n = 6 and n = 5. 
Equipped with this experience, we let Rate work out the fourth coordinate: 

In[27]:= Rate[0,0,0,-1,-6,-20] 


2 

-((-3 + iO) (-2 + iO) (-1 + iO)) 

Out [27]= {-} 

12 

After replacement of iO by n — 5, this is —{n — 8){n — 7)^(n — 6)/12. Let us summarise 
our results so far: the hrst hve coordinates of our vector in the kernel of are 

_2 -(2n-5) (^-4)(2n-8) (n - 6 )(n - 5)(2u - 11) 

(n-8)(n-7)(u-6)(2u-14) 

12 ■ 

I would say, there is a clear pattern emerging: the s-th coordinate is equal to 
{n — 2s)s-i {2n — 3s — 2) (2n — 3s — 2) (n — 2s)s 

s! (n — s — 1 ) s! 

Denoting the above expression by f{n, s), the vector 

{f{n, 0), f{n, 1 ),..., f{n, n - 2 ), 1) 

is apparently in the kernel of for n >2. To prove this, we have to show that 


n—2 


E 


(2n — 3s — 2) (n — 2s)s 
(n — s — 1 ) (s)! 


• {-y -2n + i + 2)s {-2y - 2n + 2i - s + 3)s {y + 2s - i + l) 2 n- 2 s -2 

= {—y — 2n + i + 2)„_i {—2y — 3n + 2i + 4)„_i. 


In jl()5] it was argued that this identity follows from a certain hypergeometric iden¬ 
tity due to Singh |lfi9j . However, for just having some proof of this identity, this 
careful literature search was not necessary. In fact, nowadays, once you write down 
a binomial or hypergeometric identity, it is already proved! One simply puts the bi- 
no mial/hypergeometric sum into the Gosper-Zeilberger algorithm (see jl44[ 1193111941 
[isn]), which outputs a recurrence for it, and then the only task is to verify that the (con¬ 
jectured) right-hand side also satishes the same recurrence, and to check the identity 
for sufficiently many initial values (which one has already done anyway while producing 
the conjecture).^® 


you may have suspected, this is again a little bit oversimplified. But not much. The 
Gosper-Zeilberger algorithm applies always to hypergeometric sums, and there are only very few 
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As I mentioned earlier, actually we need more vectors in the kernel. However, this 
is not difficult. Take a closer look, and you will see that the pattern persists (set 
c[n — 1] = 0 in the vector for H[n], etc.). It will take you no time to work out a 
full-fledged conjecture for \_{n -|- 1)/3J linear independent vectors in the kernel of M„. 

I do not want to go through Steps (S2) and (S3), that is, the degree calculation and 
the computation of the constant. As it turns out, to do this conveniently you need 
to introduce more variables in the determinant in (I13D. Once you do this, everything 
works out very smoothly. I refer the reader to ma for these details. 

Now, let us come back to our determinant, the determinant of (O, and apply 
“identihcation of factors” to it.^^ To begin with, here is bad news: “identihcation 
of factors” crucially requires the existence of indeterminates. But, where are they in 
(EHD ? If we look at the definition of the matrix (Q, which, in the end, depends 
on the auxiliary functions fo,fi,go,gi defined in (USD, then we see that there are no 
indeterminates at all. Everything is (integral) numbers. So, to get even started, we 
need to introduce indeterminates in a way such that the more general determinant 
would still factor “nicely.” We do not have much guidance. Maybe, since we already 
made the abbreviation N{k) = Ak{4:k -l- 1), we should replace N{k) by X? Okay, let us 
try this, that is, let us put 

/o(J)=J(-4)^ 

/l(J) = -(4J + 2)(-4)^ 

go{j) = {X -j), 

g,{j) = -{AX-Aj-2) (4.5) 

instead of (ESD- Let us compute the new determinant for k = 2. We program the new 
functions /o,/i, 5'o, 


Hi 

O 

1 — 1 

t- , 

j-] 

:= j^C-dj^k; 


t_ , 

j-] 

:= -(2 + 4*j)*(-4)^k; 

1_1 

o 

hO 

t_, 

j-] 

:= (X - j); 


t-, 

j-] 

:= (-4*X + 2 + 4*j) 


we enter the new determinant for A: = 2, 
In [29]:= Factor[Det[A[2]]] 


binomial sums where it does not apply. (For the sake of completeness, I mention that there are 
also several algorithms available to deal with multi-sums, see jsiiiini. These do, however, rather 
quickly challenge the resources of today’s computers.) Maple implementations written by Doron 
Zeilberger are available from http://www.math.temple.edu/~zeilberg, those written by Frederic 
Chyzak are available from http://algo.inria.fr/chyzak/mgfun.html, Mathematica implementa¬ 
tions written by Peter Paule, Axel Riese, Markus Schorn, Kurt Wegschaider are available from 
http://www.rise.uni-linz.ac.at/research/combinat/risc/software 

^“^What I describe in the sequel is, except for very few excursions that ended up in a dead end, and 
which are therefore omitted here, the way how the determinant evaluation was found. 
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and, after a waiting time of more than 15 minntes,^® we obtain 

Out[29]= -1406399608474882323154910525986578515918369681041517636\ 

2 

> 11783762359972003840000000 (-64 + X) (-48 + X) (-40 + X) 

3 4 5 6 7 

> (-32 + X) (-24 + X) (-16 + X) (-8 + X) X 


2 


> (9653078694297600 - 916000657637376 X + 36130368757760 X 


3 4 5 6 

> - 758218948608 X + 8928558848 X - 55938432 X + 145673 X ) 

Not bad. There are many factors which are linear in X. (This is what we were after.) 
However, the irredncible polynomial of degree 6 gives ns some headache. (The degrees 
of the irredncible part of the polynomial will grow qnickly with k.) How are we going 
to guess what this factor could be, and, even more daunting, even if we should be able 
to come up with a guess, how would we go about to prove it? 

So, maybe we should modify our choice of how to introduce indeterminates into the 
matrix. In fact, we overlooked something: maybe, in a hidden manner, the variable X is 
also there at other places in (ESD, that is, when X is specialised to N{k) = Ak{Ak + 1) 
at these places it becomes invisible. More specihcally, maybe we should insert the 
difference X — Ak{Ak + 1) in the dehnitions of /o and /i (which would disappear for 
X = Ak{Ak + 1)). So, maybe we should try: 

fo{j) = {Ak{Ak + 1) - X + j){-A)\ 

Mj) = -{lQk{Ak + 1) - 4X + 2 + 4J)(-4)^ 

9o{j) = {X -j), 
g^{j) = -iAX-Aj-2), 

Okay, let us modify our computer program accordingly. 


^®which I use to explain why our computer needs so long to calculate this determinant of a very 
sparse matrix of size 16 • 2^ = 64: isn’t it true that, nowadays, determinants of matrices with several 
hundreds of rows and columns can be calculated without the slightest difficulty (particularly if they are 
very sparse)? Well, we should not forget that this is true for determinants of matrices with numerical 
entries. However, our matrix (13.811 with the modified definitions (14.511 of /o,/i,5oil?i has now entries 
which are polynomials in X. Hence, when our computer algebra program applies (internally) some 
elimination algorithm to compute the determinant, huge rational expressions will slowly build up and 
will slow down the calculations (and, at times, will make our computer crash ...). As I learn from 
Dave Saunders, Maple and Mathematica do currently in fact not use the best known algorithms for 
dealing with determinants of matrices with polynomial entries. (This may have to do with the fact 
that the developers try to optimise the algorithms for numerical determinants in the first case.) It is 
known how to avoid the expression swell and compute polynomial matrix determinants in time about 
mn^, where n is the dimension of the matrix and m is the bit length of the determinant (roughly, in 
univariate case, m is degree times maximum coefficient length). 
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l-h 

0 

1—1 

1 

t_, 

j-] 

:= (4*k*(4*k + 1) - X + j)*(-4)^k; 

fl[k_, 

t_, 

j-] 

:= -(4*4*k*(4*k + 1) - 4*X + 2 + 4*j)*(-4)^k; 

gO[k_, 

t_, 

j-] 

:= (X - j); 

gl[k_, 

t_, 

j-] 

:= (-4*X + 2 + 4*j) 


and let us compute the new determinant for k = 2: 

In[31] := Factor [Det [A [2] ] ] 

This makes us wait for another 15 minutes, after which we are rewarded with: 

Out[31]= -296777975397624679901369809794412104454134763494070841\ 

> 1155365196124754770317472271790417634937439881166252558632\ 

> 616674197504000000000 (-141 + 2 X) (-139 + 2 X) (-137 + 2 X) 

> (-135 + 2 X) (-133 + 2 X) (-131 + 2 X) (-129 + 2 X) 

Excellent! There is no big irreducible polynomial anymore. Everything is linear factors 
in X. But, wait, there is still a problem: in the end (recall Step (S2)!) we will have 
to compare the degrees of the determinant and of the right-hand side as polynomials 
in X. If we expand the determinant according to its dehnition, then the conclusion is 
that the degree of the determinant is bounded above by 16/c^ — 1 , which, for k = 2 is 
equal to 31. The right-hand side polynomial however which we computed above has 
degree 7. This is a big gap! 

I skip some other things (ending up in dead ends ...) that we tried at this point. 
Altogether they pointed to the fact that, apparently, one indeterminate is not sufficient. 
Perhaps it is a good idea to “diversify” the variable X, that is, to make two variables, 
Xi and X 2 , out of X: 

fo{j) = m4k + l)-X2+j){-4)\ 

/i(j) = -{16k{4k -F 1 ) - 4Xi + 2 + 4J)(-4)^ 

9o{j) = (^2 - j), 
g^{j) = -{4X,-4j-2). 

We program this, 

In[32]:= f0[k_, t_, j_] := (4*k*(4*k + 1) - X[2] + j)*(-4)^k; 

fl[k_, t_, j.] := -(4*4*k*(4*k + 1) - 4*X[1] + 2 + 4*j)*(-4)^k; 

g0[k_, t_, j.] := (X[2] - j); 

gl[k_, t_, j.] := (-4*X[1] + 2 + 4*j) 

and, in order to avoid overstraining our computer, compute this time the new determi¬ 
nant for k = 1: 

In[33] := Factor [Det [A [1] ] ] 

After some minutes there appears 

Out[33]= 3242591731706757120000 (-37 + 2 X[l]) (-35 + 2 X[l]) 

3 2 

> (-33 + 2 X[l]) (1 + 2 X[l] - 2 X[2]) (3 + 2 X[l] - 2 X [2]) 
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> (5 + 2 X[l] - 2 X[2]) 

on the computer screen. On the positive side: the determinant still factors completely 
into linear factors, something which we could not expect a priori. Moreover, the degree 
(in Xi and X 2 ) has increased, it is now equal to 9 although we were only computing 
the determinant for k = 1. However, a gap remains, the degree should be 16k^ — 1 = 15 
if k = 1. 

Thus, it may be wise to introduce another genuine variable, Y. For example, we may 
think of simply homogenising the dehnitions of fo,fi,go,gi'- 

foU) = + 1)Y -X 2 + JT)(-4)^ 

/i(j) = -{16k{4k + 1)Y - 4X, + (2 + 4J)H)(-4)^ 

go{j) = {X 2 -jY), 

g,ij) = -iAX,-iAj + 2)Y). 

We program this. 


f0[k_, 

t_, 

j-] 

:= (4*k*(4*k + 1)*Y - X[2] + j*Y)*(-4)^k; 

fl[k_, 

t_, 

j-] 

:= -(4*4*k*(4*k + 1)*Y - 4*X[1] + 

(2 + 4*j)*Y)*(-4)^k; 

gO [k_, 

t_, 

j-] 

:= (X[2] - j*Y); 


t-, 

j-] 

:= (-4*X[1] + (2 + 4*j)*Y) 


In[35] := Factor [Det [A [1] ] ] 

we wait for some more minutes, and we obtain 

6 

Out[35]= -3242591731706757120000 Y (33 Y - 2 X [1]) (35 Y - 2 X[l]) 

3 2 

> (37 Y - 2 X[l]) (Y + 2 X[l] - 2 X [2] ) (3 Y + 2 X[l] - 2 X [2] ) 

> (5 Y + 2 X[l] - 2 X[2]) 

Great! The degree in Xi,X 2 ,Y is 15, as it should be! 

At this point, one becomes greedy. The more variables we have, the easier will be 
the proof. We “diversify” the variables Xi,X 2 ,Y, that is, we make them Xi^t, X 2 ^t,Yt 
if they appear in the blocks Ft or Gt, respectively, t = 1,2,... ,4/c (cf. fld.Hjl and the 
Mathematica code for the precise meaning of this dehnition): 

/o(j) = iAk{Ak + l)Yt-X 2 ,t + jYt){-A)\ 

/i(j) = -{16k{4k + l)y; - 4W,i + (2 + 4j)Yt){-4f, 

go{j) = {X2,t-jYt), 

gtij) = -i4Xt,t-i4j + 2)Yt). (4.6) 

Now there are so many variables so that there is no way to do the factorisation of 
the new determinant for fc = 1 on the computer unless one plays tricks (which we 
did).^® But let us pretend that we are able to do it: 


^®See Footnote HU for the explanation of the complexity problem. “Playing tricks” would mean to 
compute the determinant for various special choices of the variables W.t, N 2 ,t, and then reconstruct 
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In [36] := f0[k_, 

, t-. 

j-] 

:= (4*k*(4*k + l)*Y[t] - X[2, t] 

+ j*Y[t])*(-4)^k; 

fl[k_, 

, t_. 

j-] 

:= -(4*4*k (4*k + l)*Y[t] - 4*X[1, t] + 
(2 + 4*j)*Y[t] )*(-4)'^k; 

gO [k_, 

, t_. 

j-] 

:= (X[2, t] - j*Y[t]); 


, t_. 

j-] 

:= (-4*X[1, t] + (2 + 4*j)*Y[t]) 


In[37]:= Factor[Det[A[1]]] 

Out[37]= 3242591731706757120000 (2 X[l, 1] - 33 Y[l]) Y[l] 

> (2 X[l, 1] - 2 X[2, 1] + Y[l]) (2 X[l, 2] - 35 Y[2]) Y[2] 

> (2 X[l, 2] - 2 X[2, 2] + Y[2]) (-2 X[2, 2] Y[l] + 2 X[l, 1] 

> Y[2] + 3 Y[l] Y[2]) (2 X[l, 3] - 37 Y [3]) Y[3] (2 X[l, 3] - 

> 2 X[2, 3] + Y[3]) (-2 X[2, 3] Y[l] + 2 X[l, 1] Y[3] + 

> 5 Y[l] Y[3]) (-2 X[2, 3] Y[2] + 2 X[l, 2] Y[3] +3 Y [2] Y[3]) 


By staring a little bit at this result (and the one that we computed for /c = 2), we 
extracted that, apparently, we have 

Ak 

det = (_l)A:-l42fc(4fc2+7fc+2)^2/c(4fc+l) ^ 

i=l 

4fc-l 

X J] (2Xi,, - (32fc2 + 2a - l)i;) 

a=l 

X H (2X2,bi;-2Xi,,n-(26-2a+l)y;n), (4.7) 

l<a<b<Ak-l 

where denotes the new general matrix given through (UHl) and (HID, and where, 
as before, {a)k is the standard notation for shifted factorials (Pochhammer symbols) 

the general result by interpolation. This is possible because we know an a priori degree bound (namely 
15) for the polynomial. However, this would become infeasible for k = 3, for example. “Playing tricks” 
then would mean to be content with an “almost sure” guess, the latter being based on features of the 
(unknown) general result that are already visible in the earlier results, and on calculations done for 
special values of the variables. For example, if we encounter determinants detM(fc), where the M(fc)’s 
are some square matrices, fc = 1,2,..., and the results for fc = 1, 2,..., fco — 1 show that x — y must 
be a factor of det M{k) to some power, then one would specialise y to some value that would make 
X — y distinct from any other linear factors containing x, and, supposing that y = 17 is such a choice, 
compute det M (/cg) with y = 17. The exact power oi x — y in the unspecialised determinant det M (fcg) 
can then be read off from the exponent of cc — 17 in the specialised one. If it should happen that it is 
also infeasible to calculate detM(fco) with x still unspecialised, then there is still a way out. In that 
case, one specialises y and x, in such a way that x — y would be a prime p that one expects not to 
occur as a prime factor in any other factor of the determinant detM(fco). The exact power ot x — y 
in the unspecialised determinant detM(fco) can then be read off from the exponent of p in the prime 
factorisation of the specialised determinant. See Subsections 15 . 71 and 15.XI and in particular Footnote 1771 
for further instances where this trick was applied. 
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explained in the statement of Theorem [Tj The special case that we need in the end to 
prove our Theorem |S1 is Xi^t = X 2 ^t = N{k) and Yt = 1. 

Now we are in business. Here is the Sketch of the proof of dUZl): 

Re (SI): For each factor of the (conjectured) result (14.71) . we hnd a linear combination 

of the rows which vanishes if the factor vanishes. (In other terms: if the indeterminates 

in the matrix are specialised so that a particular factor vanishes, we hnd a vector in the 
kernel of the transpose of the specialised matrix.) For example, to explain the factor 
(2Xi^i — (32fc^ + l)Fi), we found: 

If Xi,i = then 


2(X2,4fc-l - (Nik) - l)y4fc-i) 

-4)H4fc+i)+i(i6A:2 + 1) + 1) 

4k 4k-2 / , , s-1 


(row 0 of H^) 




1 + Uk 2Xi,£ 


s=0 t=0 


4.fc iii6p + i_4£;t 11 

l=Q e=4k-t 


(32fc2 + 2i- l)Yt A 

(16A;2 + £ - l)Yt-i J 


■ (row (16k‘^ — (4k — l)s — t — 1) of A^) = 0, (4.8) 


as is easy to verify. (Since the coefficients of the various rows in fj4.8|l are rational 
functions in the indeterminates Xi^t, X 2 ^t^Yt, they are rather easy to work out from 
computer data. One does not even need Rate ...) 


Re (S2): The total degree in the Xi^t’s, X 2 ,ts, Yfs of the product on the right-hand 
side of dUD is 16fc^ — 1. As we already remarked earlier, the degree of the determinant 
is at most 16fc^ — 1. Hence, the determinant is equal to the product times, possibly, a 
constant. 

Re (S3): For the evaluation of the constant, we compare coefficients of 

Y‘ik Y‘^k-1 -^2 

-'2 ■ ■ ■ ■'4fc-l ■ 

After some rehection, it turns out that the constant is equal to a determinant of the 
same form, that is, of the form (EHD, but with auxiliary functions 

f,(j) = (N(k)+j)(-4)\ 

^(J)=4(-4)^ 

9o(j) = -J, 

9i(j) = -4. (4.9) 


What a set-back! It seems that we are in the same situation as at the very beginning. 
We started with the determinant of the matrix (EHD with auxiliary functions (El), and 
we ended up with the same type of determinant, with auxiliary functions (El. There 
is little hope though: the functions in (El are somewhat simpler as those in (El- 
Nevertheless, we have to play the same game again; that is, if we want to apply the 
method of identihcation of factors, then we have to introduce indeterminates. Skipping 







28 


C. KRATTENTHALER 


the experimental part, we came up with 

/o(J) = (^* + J)(-4)^ 

/i(j) = 4(-4)^X„ 

^o(j) = -J, 

9iij) = -4Xt, 

where t has the same meaning as before in (USD- 
computer calculations suggested that apparently 

4fe-l 

detAl^ = + 1)! JJ 

a=l 

The special case that we need is Zt = N{k) and Xt 

So, we apply again the method of identihcation of factors. Everything runs smoothly 
(except that the details of the verihcation of the factors are somewhat more unpleasant 
here). When we come hnally to the point that we want to determine the constant, it 
turns out that the constant is equal to — no surprise anymore — the determinant of a 
matrix of the form (O with auxiliary functions 

/o(j) = (-4)^ 

/l(J)=4(-4)^ 

^o(j) = 0, 

^i(j) = -4. 

Now, is this good or bad news? In other words, while painfully working through the 
steps of “identihcation of factors,” will we forever continue producing new determinants 
of the form which we must again handle by the same method? To give it away: 

this is indeed very good news. The function go{j) vanishes identically! It makes it 
possible that now Method 0 (= do some row and column manipulations) works. (See 
13 for the details.) We are — hnally — done with the proof of (HTj) . and, since the 
right-hand side does not vanish for Xi^t = = N{k) and Y) = 1, with the proof of 

Theorem IHl ,—. 


Denoting the new matrix by , 


a—1 


^4fc+i-an(z,-46A:) ). (4.10) 

= 1 . 


b=0 


5. More determinant evaluations 

This section complements the list of known determinant evaluations given in Section 3 
of |l()0j . I list here several determinant evaluations which I believe are interesting or 
attractive (and, in the ideal case, both), that have appeared since jlOOj . or that I failed 
to mention in jl09j . I also include several conjectures and open problems, some of 
them old, some of them new. As in uni, each evaluation is accompanied by some 
remarks providing information on the context in which it arose. Again, the selection 
of determinant evaluations presented reflects totally my taste, which must be blamed 
in the case of any shortcomings. The order of presentation follows loosely the order of 
presentation of determinants in uni- 
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5.1. More basic determinant evaluations. I begin with two determinant eval¬ 
uations belonging to the category “standard determinants” (see Section 2.1 in unni)- 
They are among those which I missed to state in uni. The reminder for inclusion here 
is the paper [i. There, Amdeberhan and Zeilberger propose an automated approach 
towards determinant evaluations via the condensation method (see “Method 2” in Sec¬ 
tion iD- They provide a list of examples which can be obtained in that way. As they 
remark at the end of the paper, all of these are special cases of Lemma 5 in uni, with 
the exception of three, namely Eqs. (8)-(10) in [i. In their turn, two of them, namely 
(8) and (9), are special cases of the following evaluation. (For (10), see Lemma ITTI 
below.) 

Lemma 9. Let P{Z) be a polynomial in Z of degree n — 1 with leading coefficient L. 
Then 

_deq(P(X. + F,)) = L“n(""b n (5.1) 

□ 


This lemma is easily proved along the lines of the standard proof of the Vandermonde 
determinant evaluation which we recalled in Section |3] (see the proof of (D) or by 
condensation. 

A multiplicative version of Lemma mi is the following. 


Lemma 10. Let P{Z) = pn-iZ^ ^ + pn-2Z^ ^ + ■ ■ ■ + Po- Then 

n—1 

j_eyyp(.xy))) = nR If (x,-x,m-Y,). 

i=0 


(5.2) 

□ 


On the other hand, identity (10) from ji can be generalised to the following Cauchy- 
type determinant evaluation. As all the identities from P, it can also be proved by the 
condensation method. 


Lemma 11. Let Oq, Oi,..., a„_i, Cq, Ci,..., c„_i, b, x and y be indeterminates. Then, 
for any positive integer n, there holds 


det 


{x + at + Cj){y + bi + cj) 


\ {x + Oi + bi + Cj) 

= 6”“^ (n - 1)! I {n-l)x + y + ^ai + ^Ci 


n—1 n—1 


2=1 2=0 


n—1 


X 


n (cj - Ci) Yliy - X - tti) JJ {{j-i)b-ai + 
0<i<j<n—l 2=1 


n—1n—1 

nn {x + tti + bi + Cj) 

i=l j=0 


(5.3) 


□ 




30 


C. KRATTENTHALER 


Speaking of Cauchy-type determinant evaluations, this brings us to a whole family 
of such evaluations which were instrumental in Kuperberg’s recent advance lira on 
the enumeration of (symmetry classes of) alternating sign matrices. The reason that 
determinants, and also Pfaffians, play an important role in this context is due to Propp’s 
discovery (described for the hrst time in Sec. 7] and exploited in |llbl I117] i that 
alternating sign matrices are in bijection with configurations in the six vertex model, and 
due to determinant and Pfaffian formulae due to Izergin [SH] and Kuperberg jll7j for 
certain multivariable partition functions of the six vertex model under various boundary 
conditions. In many cases, this leads to determinants which are, or are similar to, 
Cauchy’s evaluation of the double alternant (see jldfil vol. Ill, p. 311] and (15.5|) below) 
or Schur’s Pfaffian version |166L pp. 226/227] of it (see (I5.7|l below). 

Let me recall that the Pfaffian Pf(2l) of a skew-symmetric (2n) x (2n) matrix A is 
dehned by 

pf(A)=n (5-4) 

TT (ij)€7r 

where the sum is over all perfect matchings tt of the complete graph on 2n vertices, 
where c(7r) is the crossing number oi tt, and where the product is over all edges {ij), 
i < j, in the matching tt (see e.g. da Sec. 2]). What links Pfaffians so closely to 
determinants is (aside from similarity of dehnitions) the fact that the Pfaffian of a 
skew-symmetric matrix is, up to sign, the square root of its determinant. That is, 
det(74) = Pf(74)^ for any skew-symmetric (2n) x (2n) matrix A (cf. jl79[ Prop. 2.2]). 
See the corresponding remarks and additional references in jl()91 Sec. 2.8]. 

The following three theorems present the relevant evaluations. They are Theo¬ 
rems 15-17 from uni. All of them are proved using identihcation of factors (see 
“Method 3” in Section 0}. The results in Theorem ini contain whole sets of indeter- 
minates, whereas the results in Theorems M and HU only have two indeterminates p 
and q, respectively three indeterminates p, q and r, in them. Identity (lEHl) is originally 
due to Laksov, Lascoux and Thorup mni and Stembridge da. independently. The 
reader must be warned that the statements in urn Theorems 15-17] are often blurred 
by typos. 


Theorem 12. Let xi,X 2 , ■ ■ ■ and yi,y 2 , • • • be indeterminates. Then, for any positive 
integer n, there hold 


det ( —i ] 

l<i,j<n \Xi + PjJ 


JJ {Xi - Xj){yi - yj) 

l< 2 <j<n 


n 


(5.5) 




JJ (1 - XiXj){l - Piyfiixj - Xi){yj - yf) 

n {xi + yfiil + Xiyfi 

n 

X n(l-2/*), (5-6) 

i=l 


det 


Xi + yj 1 Xiyj 
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Pf f- - 

l<i,j<2n \Xi + Xj J 


= n 

l<i<j<2n 


Ob 2 Ob j 

Xi + Xj' 


(5.7) 


Pf 


l<i,j'< 2 n 


/ Xj - Xj \ 
\1- XiXjJ 


n 

l<i<j<2n 


Ob 2 Ob J 

1 — XiXj 


(5.8) 

□ 


Theorem 13. Let p and q he indeterminates. Then, for any positive integer n, there 
hold 


det 

l<2j<n 


q-n+j-i _ q-{n+j-i) 
pn+j-i _ p-in+j-i) 


ijp — p {qpd ^ — q ^p 

l<2,7<n 

(5.9) 


det 


1<^ j<n 


/ qj-i ^ q-U-i)\ 

Vp^-* + p-b-b J 


2 ” n (r-'-p-’*-) n ( 9 p^" - tv-'-"') 

(_1)(S_!!L!__ 

n (p’-+p-'+‘) 

(5.10) 


det 


qU+j+i _ q-{n+j+i) qTi+j-i _ q-{n+j-i) 


l<i,i<n ypn+j+i — p-{n+j+i) pn+j-i _ p-in+j-i) 

Ul<i<j<2niP’~' - ni<Li<2n+l(gp^“* - 

■ 2|i 


n (p 

l<2,ji<n 


n+j-i _ ^-{n+j-p-^^pn+j+i _ ^-{n+j+i) 


P 


(5.11) 


p 


qj+i + q-U+i) qj-i + q-i+i 


^<i,j<n + p~b+b pi~^ + p-ij-b 

2”- 


p 


n ^P 


l<2<j'<n 


= -1 u 


l<iJ<2n-\-l 

2 tL 2 |i 


JJ (P^-* + p-(l-*))(p^'+* + 

l<2,ji<n 


(5.12) 


□ 



32 


C. KRATTENTHALER 


Theorem 14. Let p, q, and r he indeterminates. Then, for any positive integer n, 
there hold 


l<iJ<2n Y (pJ”* — ^“(-5“*)) / 

Yl {p^+^-^ - p-^^+^-^i) 

l< 2 ,ji<n 


(5.13) 


/ ni+i — n-U+i) nJ-i — Q-lj-i) 

pf I ('pj+* _ ( - _ - _ - _ - _ 

idk^nV^ ^ ^ \pi+i - p-ii+i) pj-i-p-U-i) 


'■j+i _ -f-U+i) j-j-i _ f-U-i) 


pj+i — p-{j+i) pj-i — p-(j-i) 


JJ (pj' 0 *)) JJ (gp, 

l<2<_7<2n l<2j'<2n+l 

21 j 




q p 


n 

l<i<j<2n 


/■+* _ p-U+i)'' 


(5.14) 


□ 


Subsequent to Kuperberg’s work, Okada una related Kuperberg’s determinants and 
Pfafiians to characters of classical groups, by coming up with rather complex, but still 
beautiful determinant identities. In particular, this allowed him to settle one more of 
the conjectured enumeration formulae on symmetry classes of alternating sign matri¬ 
ces. Generalising even further, Ishikawa, Okada, Tagawa and Zeng jS3] have found more 
such determinant identities. Putting them into the framework of certain special rep¬ 
resentations of the symmetric group, Lascoux jl22j has clarified the mechanism which 
gives rise to these identities. 

The next six determinant lemmas are corollaries of elliptic determinant evaluations 
due to Rosengren and Schlosser |15t)j . (The latter will be addressed later in Subsec¬ 
tion Em) They partly extend the fundamental determinant lemmas in jl()91 Sec. 2.2]. 
For the statements of the lemmas, we need the notion of a norm of a polynomial 
oq + aiz + ■ ■ ■ + , which we define to be the reciprocal of the product of its roots, 

or, more explicitly, as (—1)*’ 0 ^/ 00 . 

If we specialise p = 0 in Lemma IHHl (15.12211 . then we obtain a determinant identity 
which generalises at the same time the Vandermonde determinant evaluation. Lemma IHl 
and Lemma cni 

Lemma 15. Let Pi, P 2 ,... ,Pn he polynomials of degree n and norm t, given by 

n—1 

k={) 
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Then 


det iPAxi)) 


(1 - tXi ■■■Xn) 


n 

l<2<J'<n 


det (aji). 

l<i<n 

0 <j<n-l 


(5.16) 

□ 


Further determinant identities which generalise other Weyl denominator formulae (cf. 
[mm Lemma 2]) could be obtained from the special case p = 0 of the other determinant 
evaluations in Lemma IHHl 

A generalisation of Lemma 6 from uni in the same spirit can be obtained by setting 
p = 0 in Theorem ESI It is given as Corollary 5.1 in |156j . 

Lemma 16. Let xi,..., Xn, oi,..., On, and t be indeterminates. For each j = 1,..., n, 
let Pj be a polynomial of degree j and norm tai ■ ■ - aj. Then there holds 


^det^ (l-akXi) 

' k=j+i 


1 — tai ■ ■ ■ OnXi ■ ■ - Xr, 

1 -t 


JJPi(l/ai) JJ aj{xj-Xi). (5.16) 

l<2<j<n 


2 = 1 


□ 


We continue with a consequence of Theorem [7^1 (see Corollary 5.3 in |156j i. The 
special case Pj-i{x) = 1, j = 1,... ,n, is Lemma A.l of |162j . which was needed in 
order to obtain an matrix inversion that played a crucial role in the derivation of 
multiple basic hypergeometric series identities. A slight generalisation was given in [1641 
Lemma A.l]. 

Lemma 17. Let xi, ... ,Xn and b be indeterminates. For each j = 1,... ,n, let Pj_i(x) 
be a polynomial in x of degree at most j — I with constant term 1, and let Q(x) = 
(1 — pix) • • ■ (1 — p„+ix). Then there holds 


Q{b) 


det 

l<2j<n 


X, 


n+l-j 


Pj-liXi) 




Q^\ 

Qib)J 


n 

= {I-hxi-■ ■ Xnyi-■ ■yn+i)W{xi-b) {Xi-Xj). (5.17) 

2 = 1 1<^<J'<72 


□ 

Pairing the (f, j)-entry in the determinant in (I5.16|l with itself, but with x* replaced 
by 1/xj, one can construct another determinant which evaluates in closed form. The 
result given below is Corollary 5.5 in |156j . It is the special case p = 0 of Theorem 1771 

Lemma 18. Let xi,..., x„, Oi,..., and Ci,..., Cn +2 be indeterminates. For each 
j = 1,... ,n, let Pj be a polynomial of degree j with norm (ci • • • c„+ 2 aj+i ■ • • a„)“h 
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Then there holds 


det 


n+2 


X, 


JJ (1 - CfcXi) Pj(Xi) JJ (l-ttkXi) 


k=l 


k=j+l 


n+2 


-X, 


n+1 


jj(i - ckXi Pj{Xi n 

k=j+l 


-1^ 


k=l 


ai-■ ■ a„ 


n^ 4 l/a.) 


X 


Xi ' ' ‘ Xfi (1 Cl * * * Cji-\-20>l ' ' ' Oju) 

n 

n n ajixi-xj){i-i/xixj). 


l<i<j<n-\-2 


i=l 




(5.18) 


□ 

It is worthwhile to state the limit case c „+2 —>■ oo of this lemma separately, in which 
case the norm constraint on the polynomials Pj drops out, but, in return, the degree of 
Pj gets lowered by one (see |1561 Cor. 5.8]). 

Lemma 19. Let xi,..., 02 ,..., an, and ci,..., Cn+i be indeterminates. For each 

j = 1, ... ,n, let Pj-i he a polynomial of degree at most j — 1. Then there holds 


n+1 n 

det xf'^ TT (1 “ ^kXi) Pj-i{xi) TT (1 - ttkXi) 

k=l k=j-\-l 


n+1 




k=l 


k=j+l 


Y[Pi-i{l/ai) JJ ( 1 -CiCj) 


2 = 1 


X 


Y[xi ^(1 - x^i) JJ aj{xi - Xj){l - 1 /xiXj). (5.19) 

l< 2 <j<n 


2=1 


□ 


Dividing both sides of (I5.19j] by 11^=2 then letting tend to cxo, i = 

2,3,...,n, we arrive at the determinant evaluation below (see |15t)l Cor. 5.11]). Its 
special case Pj-i{x) = 1, j = 1,..., n, is Lemma A.11 of |162j . needed there in order to 
obtain a Cn matrix inversion, which was later applied in ng to derive multiple q-Abel 
and q-Rothe summations. 
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Lemma 20. Let Ti, ..., Xn, and ci,..., Cn+i be indeterminates. For each j = 1,... ,n, 
let Pj-i be a polynomial of degree at most j — 1. Then there holds 


n+1 

j ]^(1 - CkXi) Pj-i[Xi) - x\ 


n+l 

- CkXf^) Pj-l{,X-^) 

k=l 
n 

Wxf^{l-xl) JJ {xj-Xi){l-l/xiXj). (5.20) 


k=l 


n (1 

l<i<_ 7 <n+l 




2=1 


2 = 1 


l<2<j<r2 


□ 


It is an attractive feature of this determinant identity that it contains, at the same 
time, the Weyl denominator formulae for the classical root systems Bn, Cn and Dn as 
special cases (cf. jlOOl Lemma 2]). This is seen by setting Pj{x) = 1 for all j, Ci = C 2 = 
■ ■ ■ = Cn-i = 0, and then c„ = 0, c„+i = —1 for the type Bn case, c„ = c„+i = 0 for the 
type Cn case, and c„ = 1, Cn+i = — 1 for the type Dn case, respectively. 

A determinant which is of completely different type, but which also belongs to the 
category of basic determinant evaluations, is the determinant of a matrix where only 
two (circular) diagonals are hlled with non-zero elements. It was applied with advantage 
in IHHl to evaluate Seott-type permanents. 

Lemma 21. Let n and r be positive integers, r < n, and Xi,X 2 , ■ ■ ■ ,Xn, yi,y 2 , ■ ■ ■ ,yn 
be indeterminates. Then, with d = gcd(r, n), we have 


/xi 0 .. 

0 X2 0 


det 


0 

yi 0 
0 1/2 


0 


0 

V 0 


0 yn—r+1 0 

0 yn—r+2 0 

0 

0 Un 


0 0 

yn—r 0 0 Xn J 

d , nid 

= Xi+(j-i)d 

i=l ^j=l 



(5.21) 


{I.e., in the matrix there are only nonzero entries along two diagonals, one of which is 
a broken diagonal.) □ 

A further basic determinant evaluation which I missed to state in [TTIT1] is the eval¬ 
uation of the determinant of a skew circulant matrix attributed to Scott cna in nssi 
p. 356]. It was in fact recently used by Fulmek in [TU] to find a closed form formula for 
the number of non-intersecting lattice paths with equally spaced starting and end points 
living on a cylinder, improving on earlier results by Forrester IS7I on the vicious walker 
model in statistical mechanics, see HOI Lemma 9]. 
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Theorem 22. Let n by a fixed positive integer, and let Oq, Oi,..., a„_i be indetermi- 


nates. 

Then 






( ao 

ai 

02 

On-2 



^n—1 

Oq 

QjI . . 

^n—3 

On-2 

det 

~0'n-2 

^n—1 

Oq 

(ln—4 

^n—3 


\ -ai 

—02 

-03 • ■ 

On—l 

Oo / 


n—1 


JJ(ao + a;2*+iai + + • ■ • + (522) 


i=0 


where u is a primitive (2n )-th root of unity. 


□ 


5.2. More coneluent determinants. Here I continue the discussion from the be¬ 
ginning of Section 3 in unni Theorems 20-24], There I presented determinant evalu¬ 
ations of matrices which, essentially, consist of several vertical strips, each of which is 
formed by taking a certain column vector and gluing it together with its derivative, its 
second derivative, etc., respectively by a similar construction where the derivative is 
replaced by a difference or g-difference operator. 

Since most of this subsection will be under the influence of the so-called “q-disease” 
we shall need the standard g-notations (a; g)^, denoting the q-shifted factorial and being 
given by (a; g)o := 1 and 

(a;g)fc := (1 - a)(l - ag) ■ ■ ■(! - aq^~^) 

if A; is a positive integer, as well as [^j^, denoting the q-binomial coefficient and being 
dehned by = 0 if fc < 0, [o]^ = 1, and 

al (l-g“)(l-g^-^)---(l-g°-^+^) 

if /c is a positive integer. Clearly we have limg^i [t]q= (^). 

^*^The distinctive symptom of this disease is to invariably raise the question “Is there also a q- 
analogue?” My epidemiological research on MathSciNet revealed that, while basically non-existent 
during the 1970s, this disease slowly spread during the 1980s, and then had a sharp increase around 
1990, when it jumped from about 20 papers per year published with the word “g-analog*” in it to over 
80 in 1995, and since then it has been roughly stable at 60-70 papers per year. In its simplest form, 
somebody who is infected by this disease takes a combinatorial identity, and replaces every occurrence 
of a positive integer n by its “q-analogue” 1 + q -\- q^ + ■ ■ ■ -\- < 7 "“^, inserts some powers of q at the 
appropriate places, and hopes that the result of these manipulations would be again an identity, thus 
constituting a “g-analogue” of the original equation. I refer the reader to the bible m for a rich source 
of g-identities, and for the right way to look at (most) combinatorial g-identities. In another form, 
given a certain set of objects of which one knows the exact number, one defines a statistics stat on 
these objects and now tries to evaluate object g®**^*^*^^ For a very instructive text following these 

lines see ESI, with emphasis on the objects being permutations. There is also an important third form 
of the disease in which one works in the ring of polynomials in variables x,y,... with coefficients being 
rational functions in g, and in which some pairs of variables satisfy commutation relations of the type 
xy = qxy. The study of such polynomial rings and algebras is often motivated by quantum groups and 
quantum algebras. The reader may want to consult Ena to learn more about this direction. While 
my description did not make this clear, the three described forms of the g-disease are indeed strongly 
inter-related. 
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The first result that I present is a g-extension of the evaluation of the confluent 
alternant due to Schendel |ltilj (cf. jl(J9l paragraph before Theorem 20]). In fact, 
Theorem 23 of jlU9j already provided a g-extension of (a generalisation of) Schendel’s 
formula. However, in [HSl Theorem 1], Johnson found a different g-extension. The 
theorem below is a slight generalisation of it. (The theorem below reduces to Johnson’s 
theorem if one puts C* = 0. For g = 1, the theorem below and |1()91 Theorem 23] 
become equivalent. To go from one determinant to the other in this special case, one 
would have to take a certain factor out of each column.) 

Theorem 23. Let n he a non-negative integer, and let A^iX) denote the nxm matrix 

(A'wUJ 

/ 0<i<n-l,0<j<m-l 

Given a composition of n, n = mi + ■ ■ ■ + mi, there holds 


C + i 
i-j 


det (^AraflXi) Am 2 {X 2 ) ■ ■ ■ AmflXi)^ 




n nn 

l<i<j<l g=l h=l 


fr TT - q^+s+h-i+ziAmr) 




. (5.23) 
□ 


In j93l Theorem 2], Johnson provides as well a confluent g-extension of the evaluation 
of Cauchy’s double alternant (j5.5|l . Already the case g = 1 seems to not have appeared 
in the literature earlier. Here, I was not able to introduce an additional parameter (as, 
for example, the C in Theorem EBI). 

Theorem 24. Let n be a non-negative integer, and let Am{X) denote the nxm matrix 

( _i_ 

V(I'i - X){Yi - qX)(Yi - q^X) ■ ■ ■ (1 ■ - qi-^X) 

Given a composition of n, n = mi -\- ■ ■ ■ -\- mi, there holds 
det (yAraflXi) AqnflX 2 ) . . . ArriflXi)^ 

nxn 

\ / mi nij 

n («- L) n If - q’-^x,) 

l<i<j<n / \l<i<j<i g=l h=l 

n i 

J] 1[{Y, - Xfl{Y, - qXfliYi - q^Xfl ■ ■ • (F, - q^^-^Xfl 





A surprising mixture between the confluent alternant and a confluent double alternant 
appears in m Theorem A.l]. Ciucu used it there in order to prove a Goulomb gas law 
and a superposition principle for the joint correlation of certain collections of holes for 
the rhombus tiling model on the triangular lattice. (His main result is in fact based on an 
even more general, and more complex, determinant evaluation, see j39l Theorem 8.1].) 
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Theorem 25. Let Si, S 2 , ■ ■ ■, Sm > 1 and ti,t 2 , ■ ■ ■ ,tn > 1 be integers. Write S = 
T = assume S > T. Let Xi,X 2 ,... ,Xm and yi,y 2 ,... ,yn be 

indeterminates. Define N to be the S x S matrix 

N =[A B] (5.25) 

whose blocks are given by 


A = 



and 


B = 


Then we have 

detN 




nT=inU{y- 


(5.26) 


(5.27) 


(5.28) 

□ 
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This theorem generalises at the same time numerous previously obtained determinant 
evaluations. It reduces of course to Cauchy’s double alternant when m = n and Si = 
S2 = ■ ■ ■ = Sm = ti = t2 = ■ ■ ■ = tn = ^- (lu that case, the submatrix B is empty.) 
It reduces to the confluent alternant for fi = ^2 = ■ ■ • = = 0 (i.e., in the case 

where the submatrix A is empty). The case m = rn, si = S 2 = ■ ■ ■ = Sm = 1, 
ti = t 2 = ■■■ = tn = r is stated as an exercise in |ld71 Ex. 42, p. 360]. Finally, a 
mixture of the double alternant and the Vandermonde determinant appeared already 
in IHi Theorem (Cauchy+)] where it was used to evaluate Scott-type permanents. This 
mixture turns out to be the special case Si = S 2 = ■ ■ ■ = = ^2 = ■ ■ ■ = = 1 

(but not necessarily m = n) of Theorem |2S1 

If S' = T (i.e., in the case where the submatrix B is empty). Theorem 1251 provides the 
evaluation of a confluent double alternant which is different from the one in Theorem 1^ 
for C = 0 and q = f. While, for the general form of Theorem 123 I was not able to hnd 
a g-analogue, I was able to find one for this special case, that is, for the case where B 
is empty. In view of the fact that there are also g-analogues for the other extreme case 
where the submatrix A is empty (namely Theorem 1221 and |1()91 Theorem 23]), I still 
suspect that a g-analogue of the general form of Theorem 123 should exist. 

Theorem 26. Let si, S2, ■ ■ ■ ,Sm > 1 and ti,t2, ■ ■ ■ An > 1 be integers such that si + 
S2 + ■ ■ ■ + Sm = + ^2 + • ■ • + ^n- L^t Xi,x2, ■ ■ ■ ,Xm and |/i, 1/2, • • •,i/n be indcterminates. 

Let A be the matrix 


A = 


oo 

[lY 


[§]g 

[hi 

[\-^]g 



(a:, .HI )*1 


{^l,yn)'‘‘ 

(^l.ynpn 


[21 

[tl] 

[11 

[21 

[tnl 

[iJ,, 

[ijg 

i 1 Jg 

[ijg 

[ijg 

1 1 Jg 

VTW 

{^i.yi)'^ 


prnv 


(^\,yn)^n + ^ 

\ 

\ 1 

rsi+tl-21 

[si-ii 

[ *1 1 

\ Si+in~21 

Lsi-ljq 


1 Sl-1 Jg 

[si-ljg 

[si-ljg 

L Sl-l Jg 

TTWW 






foi 

IE 


[01 

[11 

rtn-ii 

loJq 

loJg 

1 0 Jg 

loJg 

loJg 

1 0 Jg 


(xm,yi)‘‘ 

)'l 

{xm,yn)^ 

{xm^yn)'"^ 

{xm.yn)^'^ 

IE 

121 

rut 

[11 

[21 

rtni 

[iJ,, 

[ijg 

1. 1 Jg 

[iJg 

[iJg 

[ 1 Jg 


(xm,yi)^ 


{xm,yny^ 

{xm,yn)'^ 


r Sm 11 

r Sm 1 

r Sm~\-tl 1 

r Sm 11 

r Sm 1 

r Sm~\~in 21 

L Sm”! J q 

L Sm — 1 J q 

L Sttj, —1 

L -5m ~1 J q 

L Sm — 1 J q 

L -Sm — 1 J q 




where {x^yY ■.= {y — x){qy — x){q^y — x) ■■■ {q^ ^y — x). Then we have 


(5.29) 
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Si Sj 


det A = ge 1) I ^Xj-q'^ ^Xi) 

\l<i<j<m g=l h=l 

m n ti Sj 


U tj 


X 


n 1111(9" "s'-nnnn 


vl<2<ji<ri. 5=1 h=l 




(5.30) 

□ 


5.3. More determinants containing derivatives and compositions oe se¬ 
ries. Inspired by formulae of Mina jl35j . Kedlaya jHlI and Strehl and Wilf jl8()j for 
determinants of matrices the entries of which being given by (coefficients of) multiple 
derivatives and compositions of formal power series (see also [TIM Lemma 16]), Chu 
embedded all these in a larger context in the remarkable systematic study (321 • He 
shows that, at the heart of these formulae, there is the Fad di Bruno formulaf^ for mul¬ 
tiple derivatives of a composition of two formal power series. Using it, he derives the 
following determinant reduction formulae (23 Theorems 4.1 and 4.2] for determinants 
of matrices containing multiple derivatives of compositions of formal power series. 

Theorem 27. Let f{x) and fkix) and Wk{x), k = 0,1,... ,n, be formal power series 
in X with coefficients in a commutative ring. Then 

f-(w,(x)il>i{f(x))^'^ = (/'(x))( “ ( ('^1' (K^O) ■ 

(5.31) 

where (f)^^\x) is short for ^(f){x). If, in addition, Wk{x) is a polynomial of degree at 
most k, k = 1,2,... ,n, then 

js<n (n ”"=("=)) .isv ('('d(/(^))) ■ 

(5.32) 

□ 

Specialising the series ffx) so that the determinants on the right-hand sides of i fTTn i 
or (133a can be evaluated, he obtains numerous nice corollaries. Possible choices are 
ffx) = exp{yix), ffx) = \og{l + yix), fi^x) = x^f or ffx) = {oi + kx)/{ct + dix). See 
(331 Cor. 4.3 and 4.4] for the corresponding results. 

Further reduction formulae and determinant evaluations from address determi¬ 
nants of matrices formed out of coefficients of iterated compositions of formal power 

^^As one can read in m, “Faa di Bruno was neither the first to state the formula that bears his 
name nor the first to prove it.” In Section 4 of that article, the author tries to trace back the roots of 
the formula. It is apparently impossible to find the author of the formula with certainty. In his book 
dSl p. 312], Arbogast describes a recursive rule how, from the top term, to generate all other terms 
in the formula. However, the explicit formula is never written down. (I am not able to verify the 
conclusions in EHl. It seems to me that the author mixes the knowledge that we have today with what 
is really written in m-) The formula appears explicitly in Lacroix’s book |TT31 p. 629], but Lacroix’s 
precise sources remain unknown. I refer the reader to El Sec. 4] and (Jl for more detailed remarks 
on the history of the formula. 


det 
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series. In order to have a convenient notation, let us write for the n-fold com¬ 

position of / with itself, 

/<">W = /(/(■ ■•(/W))). 

with n occurrences of / on the right-hand side. Chu shows [23 Sec. 1.4] that it is possible 
to extend this n-fold composition to values of n other than non-negative integers. This 
given. Theorems 4.6 and 4.7 from |23I read as follows. 

Theorem 28. Let f{x) = x + Em= 2 9{x) = J2m=i9mX"^ and Wk{x), k = 
1, 2,..., n, he formal power series with coefficients in some commutative ring. Then 

^det^^(^[x^]wj{x)f^^^\g{x))^ = f^"\[ " ^ ( 11 (%'“ 2 /^) ) ’ 

\/c=l / \l^i<j<n / 

where [x^]h{x) denotes the coefficient ofx^ in the series h{x). If in addition, Wn{Q) = 0, 
then 


det 


{[x^^^]wj{x)f^y^'>{g{x))^ 



Furthermore, we have 



(5.34) 


(5.35) 

□ 


5.4. More on Hankel determinant evaluations. Section 2.7 of |lU9j was de¬ 
voted to Hankel determinants. There, I tried to convince the reader that, whenever 
you think that a certain Hankel determinant evaluates nicely, then the explanation will 
be (sometimes more sometimes less) hidden in the theory of continued fractions and 
orthogonal polynomials. In retrospect, it seems that the success of this try was mixed. 
Since readers are always right, this has to be blamed entirely on myself, and, indeed, 
the purpose of the present subsection is to rectify some shortcomings from then. 

Roughly speaking, I explained in uni that, given a Hankel determinant 

^^.det (fii+j), (5.36) 

0<2j<n—1 

to hnd its evaluation one should expand the generating function of the sequence of 
coefficients {fik)k>o in terms of a continued fraction, respectively hnd the sequence of 
orthogonal polynomials {pn{x))n>o with moments /ifc, fc = 0,1,..., and then the value of 
the Hankel determinant ()5.36|) can be read off the coefficients of the continued fraction, 
respectively from the recursion coefficients of the orthogonal polynomials. What I 
missed to state is that the knowledge of the orthogonal polynomials makes it also 
possible to hnd the value of the Hankel determinants which start with pi and fi 2 , 
respectively (instead of po)- In the theorem below I summarise the results that were 
already discussed in jl()9j (for which classical references are |1861 Theorem 51.1] or |185l 
Cor. 6, (19), on p. IV-17; Proposition 1, (7), on p. V-5]), and I add the two missing 


ones. 
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Theorem 29. Let {fik)k>o be a sequence with generating function written 

in the form 

OO 

/^o 


Then 




1 + aox 


bix^ 


(5.37) 


1 + aix — 


b2X^ 


1 + a2X 


_..det Ui+j) = ^T^b^ b^ ■ ■ ■b^_^bn-i. 

0<2 j<n—1 


(5.38) 


Let {pn{.x))n>Q be a sequence ofmonic polynomials, the polynomialPn{x) having degree 
n, which is orthogonal with respect to some functional L, that is, L{pm{x)pn{x)) = 
dm,nCn, whcrc the Cn’s are some non-zero constants and 6m,n Is the Kronecker delta. Let 


Pn+l{x) = {an + x)pn{x) - bnPn-l{x) (5.39) 

be the corresponding three-term recurrence which is guaranteed by Favard’s theorem. 
Then the generating function fof the moments pk = L{x^) satisfies (j5.37|l 

with the ai’s and bi’s being the coefficients in the three-term recurrence (j5.39|l . In 
particular, the Hankel determinant evaluation (IE3HI) holds, with the bi’s from the three- 
term recurrence (ICTll . 

If {qn) n>o is the sequence recursively defined by qo = 1, qi = —ag, and 

Qn+l (^nQn ^n?n—1) 

then in the situation above we also have 

^ ..det (/ii+i+i) = F^bf-%-^ ■ ■ ■ bl_2bn-iqn (5.40) 

0<2 j<n—1 

and 

n 

.det (/ii+j+ 2 ) = k^obi~^b'^~^ ■ ■ ■ bl_2bn-i ^ qlbk+i ■ ■ ■ bn-ibn- (5.41) 

0<2,?<n—1 ^ 

“ “ /c=0 

□ 

I did not find a reference for (lonii and (ICTD . These two identities follow how¬ 
ever easily from Viennot’s combinatorial model jl85j for orthogonal polynomials and 
Hankel determinants of moments. More precisely, in this theory the moments pk are 
certain generating functions for Motzkin paths, and, due to Theorem IHl the Hankel 
determinants deto<ij<n-i(hi+i+m) are generating functions for families (Pi, P 2 ,..., P„) 
of non-intersecting Motzkin paths, P* running from (—i, 0) to {j +m,0). In the case 
m = 0, it is explained in jl851 Ch. IV] how to find the corresponding Hankel determi¬ 
nant evaluation ()5.38p using this combinatorial model. The idea is that in that case 
there is a unique family of non-intersecting Motzkin paths, and its weight gives the 
right-hand side of (I5.38|l . If m = 1 or m = 2 one can argue similarly. The paths are 
uniquely determined with the exception of their portions in the strip 0 < a; < m. The 
various possibilities that one has there then yield the right-hand sides of fPninil and 
(|5.4ip . 

Since there are so many explicit families of orthogonal polynomials, and, hence, so 
many ways to apply the above theorem, I listed only a few standard Hankel determinant 
evaluations explicitly in |l()9j . I did append a long list of references and sketched in 
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which ways these give rise to more Hankel determinant evaluations. Apparently, these 
remarks were at times too cryptic, in particular concerning the theme “orthogonal 
polynomials as moments.’’'’ This is treated systematically in the two papers HSIEZI by 
Ismail and Stanton. There it is shown that certain classical polynomials {rn{x))n>o, 
such as, for example, the Laguerre polynomials, the Meixner polynomials, or the Al- 
Salam-Chihara polynomials (but there are others as well, see [13 Eli), are moments 
of other families of classical orthogonal polynomials. Thus, application of Theorem 123 
with /r„ = rn{x) immediately tells that the evaluations of the corresponding Hankel 
determinants 

det {ri+j{x)) (5.42) 

(and also the higher ones in (I5.4()|l and (15.4111 i are known. In particular, the explicit 
forms can be extracted from the coefficients of the three-term recursions for these other 
families of orthogonal polynomials. Thus, whenever you encounter a determinant of 
the form ()5.42|1 . you must check whether {rn{x))n>o is a family of orthogonal polyno¬ 
mials (which, as I explained in jl()9j . one does by consulting the compendium jl()2j of 
hypergeometric orthogonal polynomials compiled by Koekoek and Swarttouw), and if 
the answer is “yes”, you will find the solution of your determinant evaluation through 
the results in |H3 EZl by applying Theorem 123 


While Theorem 123 describes in detail the connexion between Hankel determinants 
and the continued fractions of the type (ESZD, which are often called J-fractions (which 
is short for Jacobi continued fractions), I missed to tell in unni that there is also a 
close relation between Hankel determinants and so-called S-fractions (which is short for 
Stieltjes continued fractions). I try to remedy this by the theorem below (cf. for example 
Theorem 7.2], where S'-fractions are called regular C-fractions). In principle, since 
S'-fractions are special cases of J-fractions (Enii in which the coefficients a* are all zero, 
the corresponding result for the Hankel determinants is in fact implied by Theorem |23 
Nevertheless, it is useful to state it separately. I am not able to give a reference for 
(j5.4b|) . but, again, it is not too difficult to derive it from Viennot’s combinatorial model 
CHSI for orthogonal polynomials and moments that was mentioned above. 


Theorem 30. Let {p^k)k>o be a sequence with generating function J2T=oLkX^ written 
in the form 


y^^l^kx’" 

k=0 


ho 


1 + 


aix 


1 + 


a2X 

1 + --- 


(5.43) 


Then 


det (h-j+i) — ho(®l® 2 )"' ^( 03 ^ 4 )"' ^ ■ (fl2n-5®2n-4)^(fl2n-3®2n-2); (5.44) 

0<2 j<n—1 

det (/Xj+j+i) = (—l)"pQa"(a 2 a 3 )"' ^( 0405 )*^ ^■■■( 02 ^- 40241 - 3 )^( 02 ^- 202 ^- 1 ), (5.45) 

0<i,j<n—l 


and 
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det (/ij+j+ 2 ) — / 4 o 0 ^( 0203 )"' ^( 0405 )*^ ^ ■ ■ ■ (a 2 n- 4 a 2 n- 3 )^(fl 2 n- 2 a 2 n-l) 

0<i,j<n—l 

X ^ ■ ■ ■ ai„. (5.46) 

0<il — l<i 2 — 2 <---<in—n<n 

□ 


Using this theorem, Tamm |1811 Theorem 3.1] observed that from Gaufi’ continued 
fraction for the ratio of two contiguous 2 Fi-series one can deduce several interesting 
binomial Hankel determinant evaluations, some of them had already been found ear¬ 
lier by Egecioglu, Redmond and Ryavec |531 Theorem 4] while working on polynomial 
Riemann hypotheses. Gessel and Xin 173 undertook a systematic analysis of this ap¬ 
proach, and they arrived at a set of 18 Hankel determinant evaluations, which I list 
as in the theorem below. They are preceded by the Hankel determinant 

evaluation (ICTI) . which appears only in inni Theorem 4]. 


Theorem 31. For any positive integer n, there hold 


det 

0<ij<n—1 




, f 1 -F 3j-H\ \ _-p-r (|)i (|)i (|)i (|)i /27\ 

+ + ^+J )) n (1)2, (1)2* UJ 

, f 1 -F 3j4\ \ _ -p-r (|)i (|)i (|)i (|)i /27\ 

o<„<n-i^3z + 3j+4V *+J + l y'y' II (1)2* (1)2* U; 

, / 1 -F 3j2\ \ _ -p-r (|)* (|)* (|)* (|)* /27\ 

o<M<n-iV3z + 3j + 2V ^+J + Wy' II (1)2* (1)2* W; 

, f 1 fSt + 3j -h 5\ \ ^ -p-r (|)* (|)* (|)* (|)* / 

o<M<n-iV3*+3j+5V*+j+211 (i)2*(i)2* U; 

/ 2 -F 3j-H\ \ ^ -p-r (|)* (|)* (|)* (|)* /27\ 

0<M<n-lU* + 3j + lU + J + l y'J ll (i)2*(|)2* U; 

, f 2 fSt -F 3j -h 4\ \ ^ -p-r (|)* (|)* (|)* (|)* 

o<^,J<n-ly3^+3J+Ay^+J+2)) 11 (1)2* (1)2* U; 


(5.47) 

(5.48) 

(5.49) 

(5.50) 

(5.51) 

(5.52) 

(5.53) 


det 

0<ij<n—l 


( 2 _ f3^ + 3J + 2\\ ^. (!)*(|)* (!)*(§)* 

V(3* + 3jl)(3i-h 3j-h 2) V ^j1 // (|)2* (|)2* 


(f)"' 

(5.54) 


det 

0<ij<n—l 


{3i + 3j + 4)(3i + 3j + 5) 


f3i + 3j + 5\\ 

V *+j + 2;; 




(5.55) 
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J . ( _ (9i + 9j + 5) _ /Si + 3j + 2\\ ^ -p-r (|)^ (|)i (|)i (^)ii / 

o<^,<n-l\{3^ + 3J + l){3^ + 3J + 2)\^+J + l JJ (1)2* (1)2. W; ’ 

(5.56) 


det 




(9i + 9j + 14) 


/3* + 3j+5\\ 


1 /' • a \ / • • \ 1 ' • r-v 

o<«j<n— 1 y (3i + 3j + 4)(3i + 3j + 5) \ i j 2 


' ' i=\ 


;)*(!)* (|)*(i)* 

(1)2* (1)2* 


Let ao = -2 and am = for m>l. Then 


n—1 


det (ai+j) = JJ(-2) 
2=0 


(l)*(-^)*(l)*(|)* f27 


,l)2i (|)2* 


0<i,j<n— 

Let bo = 10 and bm = om +2 C"m^) lr > 1. Then 

,(|)*(^)*(D*(t)* ^27 



2i 


det (bi+j) = TT 10- 


i=0 


J)2* (|)2i 


/^3i + 3j + 5\ \ _ JJ" 
7 ' i=0 


n /4 


(|)*(|)*(|)*(I)* ^27^ 


:i)2*(|)2* 


Furthermore, 

det , , 

\3i + 3j + 5 \ i + j + 1 

Let Co = I and Cm = form>l. Then 

det fr 



2i 


i=0 


:i)2*(|)2* 



Let do = -5 and dm = ^om+mm+ 2 ) (" form>l. Then 

n—1 

det {d,+j) = n(-5) 

0<2j<n —1 


(!)*(-!)* (§)*(!)* ^27 


i=0 


J)2. (|)2i 



2i 


Furthermore, 

f _^+ 3j + 5\ \ ^ -p-r 2 (|)* (§)* (|)i (|)* 

o<*,i<n-i \ {3i + 3j + 4) (3i + 3j + 5) \ i + j + 1 / / 


:i)2*(i)2* 


Let Co = 14 and (3-+2) for m>l. Then 

n-l /1 

det {ei+j) = Y[U 


Cs)^i|)^i|)^Ci)^ ^27 


0<i,j<n— 

Furthermore, 

( 2{9i + 9j + 14) 


i=0 


J)2i (|)2i 



2i 


det 


0<i,j<n-i y {3i + 3j + 4)(3i + 3j + 5) \ i + j + 1 / / 


/3. + 3j+ 5\\ A,(i)*(i)*(| 

I 11^ 


)^il 


i=l 


5' 

3n \o, 

2)2i I 2 ) 2 * 



(5.57) 


(5.58) 


(5.59) 

(5.60) 

(5.61) 


(5.62) 



(5.65) 











46 


C. KRATTENTHALER 


□ 


Some of the numbers appearing on the right-hand sides of the formulae in this theorem 
have combinatorial signihcance, although no intrinsic explanation is known why this is 
the case. More precisely, the numbers on the right-hand sides of ()5.48|1 . (Ib.blj) and 
()5.52|) count cyclically symmetric transpose-complementary plane partitions (cf. jld2j 
and 0). whereas those on the right-hand sides of (j5.47|l . (j5.49|l . (I5.5()|l and (15.5411 
count vertically symmetric alternating sign matrices (cf. jll7j and 0). 

In j5d[ Theorem 4], Egecioglu, Redmond and Ryavec prove also the following common 
generalisation of (lOnil and (15.5011 . (The hrst identity is the special case a; = 0, while 
the second is the special case a; = 1 of the following theorem.) 


Theorem 32. Form > 0, let Sm{x) = '^hen, for any positive 

integer n, there holds 


n—1 


det ( 

0<ij<n—l 


^i+j 


w)=n 


i=0 




(5.66) 


□ 


As I mentioned above, in [inni I only stated a few special Hankel determinant eval¬ 
uations explicitly, because there are too many ways to apply Theorems EHl and EHl I 
realise, however, that I should have stated the evaluation of the Hankel determinant of 
Catalan numbers there. I make this up now by doing this in the theorem below. I did 
not do it then because orthogonal polynomials are not needed for its evaluation (the 
orthogonal polynomials which are tied to Catalan numbers as moments are Chebyshev 
polynomials^ but, via Theorems EHl and IHUl one would only cover the cases m = 0,1,2 
in the theorem below). In fact, the Catalan number can be alternatively 

written as = (—and therefore the Hankel determinant evaluation be¬ 
low follows from una Theorem 26, (3.12)]. This latter observation shows that even a 
more general determinant, namely deto<ij<n-i .+j), can be evaluated in closed form. 
For historical remarks on this ubiquitous determinant see m paragraph before the 
Appendix]. 


Theorem 33. 


0<2 j<n—1 


n 


2n i j 
i+j 


(5.67) 

□ 


As in noni, let me conclude the part on Hankel determinants by pointing the reader 
to further papers containing interesting results on them, high-lighting sometimes the 
point of view of orthogonal polynomials that I explained above, sometimes a combi¬ 
natorial point of view. The hrst point of view is put forward in jlblj (see jllOj for 
the solution of the conjectures in that paper) in order to present Hankel determinant 
evaluations of matrices with hypergeometric 2 F 1 -senes as entries. The orthogonal poly¬ 
nomials approach is also used in m to show that a certain Hankel determinant dehned 
by Catalan numbers evaluates to Fibonacci numbers. In jT3] , one hnds Hankel determi¬ 
nant evaluations involving generalisations of the Bernoulli numbers. The combinatorial 
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point of view dominates in m ESI ESI IS3IMI, where Hankel determinants involving 
q-Catalan numbers, q-Stirling numbers, and q-Fibonacci numbers are considered. 

A very interesting new direction, which seems to have much potential, is opened up 
by Luque and Thibon in |127j . They show that Selberg-type integrals can be evalu¬ 
ated by means of Hankel hyperdeterminants, and they prove many hyperdeterminant 
generalisations of classical Hankel determinant evaluations. 

At last, (but certainly not least!), I want to draw the reader’s attention to Lascoux’s 
“unorthodox” approach to Hankel determinants and orthogonal polynomials through 
symmetric functions which he presents in detail in nza Ch. 4, 5, 8]. In particular. 
Theorem ESI Eq- (I5.38j) are the contents of Theorem 8.3.1 in imi (see also the end of 
Section 5.3 there), and Theorem Eqs. (I5.44|l and (I5.45|) are the contents of Theo¬ 
rem 4.2.1 in nm, The usefulness of this symmetric function approach is, for example, 
demonstrated in |H2] in order to evaluate Hankel determinants of matrices the entries 
of which are Rogers-Szego, respectively Meixner polynomials. 


5.5. More binomial determinants. A vast part of Section 3 in jlObj is occupied by 
binomial determinants. As I mentioned in Section|21of the present article, an extremely 
rich source for binomial determinants is rhombus tiling enumeration. I want to present 
here some which did not already appear in uni. 

To begin with, I want to remind the reader of an old problem posed by Andrews in 
im p. 105]. The determinant in this problem is a variation of a determinant which 
enumerates cyclically symmetric plane partitions and descending plane partitions, which 
was evaluated by Andrews in PI (see also uni Theorem 32]; the latter determinant 
arises from the one in ()5.08j) by replacing j -|- 1 by j in the bottom of the binomial 
coefficient). 

Problem 34. Evaluate the determinant 

T>i(n) := det t 

o<*j<n-i y y J -|- 1 J J 


where 6ij is the Kronecker delta. In particular, show that 

Di{2n) _ 2^^ (f + n) (f + 2n + 

Di{2n - 1) {n)n (-f -2n+\) |-(„_ 2 )/ 2 ] 


(5.69) 


could easily be that it is the “modern” treatment of the theory which must be labelled with 
the attribute “unorthodox.” As Lascoux documents in inn, in his treatment he follows the tradition 
of great masters such as Cauchy, Jacobi or Wrohski ... 
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The determinants -Di(n) are rather intriguing. Here are the hrst few values: 

-Di(l) = h + 1, 

-Oi(2) = (h + l)(h + 2), 

-Di(3) = + 1)(/^ + 2)(h + 3)(/4 + 14), 

-Di(4) = + l)(h + 2)(/4 + 3)(/i + 4)(/i + 9)(/i + 14), 

-Di(5) = (/i + l)(/x + 2)(/i + 3)(/i + 4)(/i + 5)(/r + 9) 

X (3432 + 722/r + 45/4^ + /i^), 

-Di(6) = g^g^gg(/i + l)(/r + 2)(/i + 3)(/i + 4)(/i + 5)(/i + 6)(/r + 8)(/r + 13) 

X (/i + 15) (3432 + 722/i + 45/4^ + /i^), 

-Di(7) = 3709X2000^^ + l)(h + 2)(/4 + 3)(/i + 4)(/x + 5)(/i + 6)(/i + 7)(/i + 8)^ 

X (/i + 13)(/x + 15)2(/i + 34)(/i3 + 47/4^ + 954/x + 5928), 

-Di(8) = 73X565080000^^ + l)(h + 2)(/4 + 3)(7t + 4)(/i + 5)(/i + 6)(/i + 7)(/r + 8)^ 

X (/i + 10)(/r + 15)^(/i + 17)(/i + 19) 

X (/i + 21) (/X + 34) (/x^ + 47/1^ + 954/x + 5928). 


“So,” these determinants factor almost completely, there is only a relatively small (in 
degree) irreducible factor which is not linear. (For example, this factor is of degree 6 
for Di{9) and Hi(10), and of degree 7 for Hi(11) and Hi(12).) Moreover, this “bigger” 
factor is always the same for Hi(2n — 1) and Hi(2rx). Not only that, the quotient which 
is predicted in fj5.09|l is at the same time a building block in the result of the evaluation 
of the determinant which enumerates the cyclically symmetric and descending plane 
partitions (see El)- All this begs for an explanation in terms of a factorisation of the 
matrix of which the determinant is taken from. In fact, for the plane partition matrix 
there is such a factorisation, due to Mills, Robbins and Rumsey una Theorem 5] (see 
also [Tuni Theorem 36]). The question is whether there is a similar one for the matrix 
in (j5.68|l . 

Inspired by this conjecture and by the variations in pni Theorems 11-13] (see {1091 
Theorem 35]) on Andrews’ original determinant evaluation in ^U], Guoce Xin (private 
communication) observed that, if we change the sign in front of the Kronecker delta in 
(EnHl), then the resulting determinant factors completely into linear factors. 

Conjecture 35. Let /x be an indeterminate and n be a non-negative integer. The 
determinant 


-6ij + 




det 

0<ij’<n—1 


(5.70) 
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is equal to 


(E\ ({n-2)/2 2 

^yn/2c2n{n+2) / A ^ 2/n/2 / 

( 2 )- 


i=0 


(2*) 


12 


X 


L(n-4)/4j / 2 

^ + 3z + - 

j=0 2/ („_4._2)/2 


n 


_fi_5! + 3. + 3 
2 2 


(n-4i-4)/2 


(5.71) 


z/n zs even, and it is equal to 


-l)(n-l)/22(«.+3)(n+l)/4 
'L(n-3)/4j 


/i — 1 


X 


n 

i=0 


(n+l)/2 

2 


{n-l)/2 

n 

i=0 


z!(* + l)! 

(22)1(2*+ 2)! 


/i 5\ ( fi 3n 3 

- + 3* + - --^-h32 

2 2/ („_4j_3)/2 V 2 2 2 y {n-4i-l)/2 


(5.72) 


z/n odd. 

In fact, it seems that also the “next” determinant, the determinant where one replaces 
j + 1 at the bottom of the binomial coefficient in (13701) by j' + 2 factors completely 
when n is odd. (It does not when n is even, though.) 

Conjecture 36. Let ja be an indeterminate. For any odd non-negative integer n there 
holds 


det 

0<ij<n—l 





= (^_\'^{n-l)/22{n-l){n+b)/A^^_^Y'^ 


f-l) 


(n+l)/2 


(2±1)! 



(n-4j-l)/2 


X 



3n 5 


(n-4*-3)/2 


(5.73) 


For the combinatorialist I add that all the determinants in Problem El and Con¬ 
jectures El and El count certain rhombus tilings, as do the original determinants in 

dnunmni. 

Alain Lascoux (private communication) did not understand why we should stop here, 
and he hinted at a parametric family of determinant evaluations into which the case of 
odd n of Conjecture El is embedded as a special case. 
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Conjecture 37. Let /x be an indeterminate. For any odd non-negative integers n and 
r there holds 


0<ij<n—l 


det I —+ 


F + i + j 
j + r 


^r-2 


(r-3)/2 


{n-2i-2)\^ 


f _ 1 \(n-r)/2i^(n^+6n-2nr+r^-4r+2)/4 j I I j | | I I _ 

^ ^ n (ri^)!2(^ + 2*)!(n + 2* + 2)! 


vi=0 


X 


X (/i — r) 
'l(n-r-2)/4j 

n 

i=0 


/i + 1 


/J —1 


(n-l)/2 


i\ (i + 1)! 


(n-r)/2 \i=i 


!!(“-’■ + ) I n (2,)! (2^ + 2)! 


f + 3* + r+- 


(n—4i—r—2)/2 


u 3n r 

- -+ - + 3Z + 1 

2 2 2 


(n—4i—r)/2 


(5.74) 


The next binomial determinant that I want to mention is, strictly speaking, not a 
determinant but a Pfaffian (see D for the dehnition). While doing {—!)-enumeration 
of sef-complementary plane partitions, Eisenkblbl m encountered an, I admit, compli¬ 
cated looking Pfaffian, 

Pf (M{mi,m2,ni,n2,a,b)), (5.75) 

where a is even and b is odd, and where 


Mij{mi,m2,ni,n2,a,h) 
L(a+6-l)/4)J 

l=l 


Hi W mi 

ip — l)/2 -I- [(i — 1)/2J —l + l) \—a/2 -f [{j — 1)/2J -|- I 


ni 


mi 


fb — l)/2 -|- [{j — 1)/2J — I + ij \—a/2 -f L(i — 1)/2J -|- I 

r(a+6-l)/4] 


1=1 


n2 


m2 


(b — l)/2 -|- [?/2j — / 1 y y— a/2 [j/2j 1 — 1 

n2 \ / m2 

(b — l)/2 -|- [j/2j — / 1 y y— a/2 [?/2j 1 — 1 


Remarkably however, experimentally this Pfaffian, Erst of all, factors completely into 
factors which are linear in the variables mi,m 2 ,ni,n 2 , but not only that, there seems 
to be complete separation, that is, each linear factor contains only one of mi, m 2 , ni,n 2 . 
One has the impression that this phenomenon should have an explanation in a factori¬ 
sation of the matrix in However, the task of Ending one does not seem to be an 

easy one in view of the “entangledness” of the parameters in the sums of the matrix 
entries. 


Problem 38.^^ Find and prove the closed form evaluation of the Pfaffian in (I77SD . 


^^Theresia Eisenkolbl has recently solved this problem in m 
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Our next determinants can be considered as shuffles of two binomial determinants. 
Let us first consider 


/ fb + c + m 


det 


-1 + ] 
h+c 
2 

^ - i + j + e 


1 < i < a 


\ 


a+l<i<a+m 


(5.76) 


In fact, if e = 0, and if a, 6, c all have the same parity, then this is exactly the determinant 
in (j2.1()|l . the evaluation of which proves Theorem |21 as we explained in Section |21 
If e = 1/2 and a has parity different from that of b and c, then the corresponding 
determinant was also evaluated in [lU], and this evaluation implied the companion result 
to TheoremElthat we mentioned immediately after the statement of the theorem. In the 
last section of [lOj, it is reported that, apparently, there are also nice closed forms for the 
determinant in (j5.76j) for e = 1 and £ = 3/2, both of which imply as well enumeration 
theorems for rhombus tilings of a hexagon with an equilateral triangle removed from its 
interior (see Conjectures 1 and 2 in m)- We reproduce the conjecture for e = 1 here, 
the one for £ = 3/2 is very similar in form. 


Conjecture 39. Let a, 6, c, m be non-negative integers, a, b, c having the same parity. 
Then for e = 1 the determinant in ()5.76|) is equal to 


lH(a 
4 


m) H(6 + m) H(c + m) H(a + 6 + c + m) 


H(a + 6 + m) H(a + c + m) H(6 + c + m) 

H(m+ [^^])H(m+ 


X 


X 


H(^ + m + 1) + m - 1) H(^ + m) 


X 


H(f + fll) H(f + fll) H(f + [|1) H(f + [|J) H(f + ) H(f + [|J) 

H(m)2 H(2±|t!!!)2 


' b+c' 


Pi{a,b,c,m), (5.77) 


H(f + [^])H(f + [^J)H(^^-1)H(^ + 1)H( 

where Pi{a,b,c,m) is the polynomial given by 

, , I (a + 6)(a + c) + 2am if a is even, 

PAa.b.c.m) = <) . 

(a + 5)(a + c) + 2(a + 6 + c + m)m if a is odd, 

and where H(n) is the hyperfactorial defined in ()2.8j) . 


Two other examples of determinants in which the upper part is given by one binomial 
matrix, while the lower part is given by a different one, arose in |43l Conjectures A.l 
and A.2]. Again, both of them seem to factor completely into linear factors, and both 
of them imply enumeration results for rhombus tilings of a certain V-shaped region. 
The right-hand sides of the (conjectured) results are the weirdest “closed” forms in 
enumeration that I am aware of.^^ We state just the hrst of the two conjectures, the 
other is very similar. 


non-trivial simplifications seem to be possible. 
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Conjecture 40. Let x, ?/, m be non-negative integers. Then the determinant 


det 


is egual to 


n 


l<i,j<m+y I I (, 


{x -\-i)\ 


) i = m -\-1,... ,m -\- y 


i+jj 


(5.78) 


m-\-y 

n 


(x + 2m — i + 1)! 


(x — i + m + 1 / + 1)! (2i — 1)! (2m + 2?/ — 2i + 1)! (m + x — ?/ + i — 1)! 


m—1 y—1 


X 2( 2 )^( 2 ) JJ d JJ i! JJ(x + 2 + |)m-2i-i - ?/ + | + 3i) j _2 

i=l j=l j>0 i>0 

X JJ(x + ^ - 1 / + ^ - 2 / + Lf J + 


i>0 


*>0 


X 


]^(x + m - |_|J + i + 1) 


i>0 


2[|J-m-27 


+[ij+i+2) 


i>0 


m-2[|J-2i-2 




JJ(x -y-h3i-h l)m+2y-4i JJ (x + m- 2/ + z + l) 3 y 


—m—4i 


X 


i=0 


i=0 


]^(x + Y ~ 2 ^)y-2i (a^ + y — 2 ”^^"^ §)y-2i-i 

y 

J^(x + i + 2)2m-2i-l 


i>0 


X 


i=0 


/ 5 ^^-• (5-79) 

[x + y + 2)m-y-l [m + X-y+ l)m+y 

Here, shifted factorials occur with positive as well as with negative indices. The con¬ 
vention with respect to which these have to he interpreted is 

{ a{a + 1) ■ • ■ (a + fc — 1) if k > 0, 

1 if k = 0, 

1/(0 — l)(a — 2 ) ■ • • (a + fc) if k < 0 . 

All products ni>o(/(*)) 9 (*) (|’^^79|1 have to interpreted as the products over all i > 0 

for which g{i) >0. 

For further conjectures of determinants of shuffles of two binomial matrices I refer the 
reader to Conjectures 1-3 in Section 4 of HU. All of them imply also enumeration results 
for rhombus tilings of hexagons. This time, these would be results about the number of 
rhombus tilings of a symmetric hexagon with some fixed rhombi on the symmetry axis. 


5.6. Determinants oe matrices with recursive entries. Binomial coefficients 
satisfy the basic recurrence of the Pascal triangle, 

Pij =Pi,j-i + Pi-ij. (5.80) 

We have seen many determinants of matrices with entries containing binomial coef¬ 
ficients in the preceding subsection and in |1()91 Sec. 3]. In jTTj, Bacher reports an 
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experimental study of determinants of matrices {Pi,j)o<i,j<n-i, where the coefficients 
Pij satisfy the recurrence (I5.8()|l (and sometimes more general recurrences), but where 
the initial conditions for pj q cind po,u * ^ 0, are different from the ones for binomial 
coefficients. He makes many interesting observations. The most intriguing one says 
that these determinants satisfy also a linear recurrence (albeit a much longer one). It 
is intriguing because it points towards the possibility of automatising determinant eval¬ 
uations^^, something that several authors (cf. e.g. [HI 11091 fl4dp have been aiming at 
(albeit, with only limited success up to now). The conjecture (and, in fact, a generalisa¬ 
tion thereof) has been proved by Petkovsek and Zakrajsek in jl92j . Still, there remains 
a large gap to fill until computers will replace humans doing determinant evaluations. 

The paper HU contained as well several pretty conjectures on closed form evaluations 
of special cases of such determinants. These were subsequently proved in jlllj . We 
state three of them in the following three theorems. The first two are proved in mu 
by working out the LU-factorisation (see “Method 1” in Section HD for the matrices of 
which the determinant is computed. The third one is derived by simple row and column 
operations. 

Theorem 41. Let {cii,j)i,j>o the sequence given by the recurrence 

ij T d' 1) hi — 

and the initial conditions = p* and = cr*, i > 0. Then 

det (ajj) = ( 1 -|-x)^ 2 ) (a;-)-p-|-O'_ (5.81) 

0<i,j<n—l ’ 

□ 

Theorem 42. Let (ajj)ij>o be the sequence given by the recurrence 

(^i,j = d ^ hi d 1, 

and the initial conditions 0*^* = 0, i > 0, Ojp = and Oq,, = — i >1. Then 

det (Oi ,) = (1 + a;)d-b^ (x + (5.82) 

0<i,i<2n-l ^ ' 

□ 

Theorem 43. Let {aij)ij>o be the sequence given by the recurrence 

ai,j = a-i-ij Uij-i + X i,j > 1, 

and the initial conditions Ojp = i and ao,i = —i, i > 0. Then 

det (oj o) = (1 -I- (5.83) 

0<i,j<2n-l V / 

□ 

Certainly, the proofs in um are not very illuminating. Neuwirth insi has looked 
more carefully into the structure of recursive sequences of the type as those in Theo¬ 
rems lUHlHl Even more generally, he looks at sequences {fi,j)i,j>o satisfying the recur¬ 
rence relation 

fi,j — d djfi,j-i d ejfi-ij-i, i,j > 1, (5.84) 

^^The reader should recall that the successful automatisation [HI iTHHi iMTinMi of the eval¬ 
uation of binomial and hypergeometric sums is fundamentally based on producing recurrences by the 
computer. 
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for some given sequences {cj)j>i, {dj)j>i, {ej)j>i. He approaches the problem by hnding 
appropriate matrix decompositions for the (inhnite) matrix {fi,j)i,j>o- In two special 
cases, he is able to apply his decomposition results to work out the LU-factorisation of 
the matrix {fi,j)i,j>o explicitly, which then yields an elegant determinant evaluation in 
both of these cases. Neuwirth’s hrst result |ld81 Theorem 5] addresses the case where 
the initial values foj satisfy a hrst order recurrence determined by the coefficients dj 
from (j5.84|l . It generalises Theorem 1^ There is no restriction on the initial values 
for i > 1. 

Theorem 44. Let {cj)j>i, (dj)j>i and {ej)j>i be given sequences, and let {fi,j)ij>o be 
the doubly indexed sequence given by the recurrence and the initial conditions 

/o,o = 1 and foj = djfoj-i, j > 1. Then 

,(Ai)= n {ei+i + Cjdi+i). (5.85) 

0 < 2 j<n—1 


Neuwirth’s second result jl88l Theorem 6] also generalises Theorem bnt in a 
different way. This time, the initial values Jqj, j > 1, are free, whereas the initial 
values fi^Q satisfy a hrst order recurrence determined by the coefficients Cj from the 
recurrence (Oil) . Below, we state its most attractive special case, in which all the c^-’s 
are identical. 

Theorem 45. Let (dj)j>i and {ej)j>i be given sequences, and let {fi,j)i,j>o be the doubly 
indexed sequence given by the recurrence fOill and the initial conditions /o,o = 1 and 
fifl = cfi-ifi, i>l. Then 

n—1 

„,.d.et {fi,j) = Y[{cdi + ei)'^-\ (5.86) 

0 < 2 ,?<n—1 

“ “ i=l 

□ 

5.7. Determinants eor signed permutations. The next class of determinants 
that we consider are determinants of matrices in which rows and columns are indexed 
by elements of reflection groups (the latter being groups generated by rehections of 
hyperplanes in real n-dimensional space; see |H1] for more information on these groups, 
and, more generally, on Coxeter groups). The prototypical example of a rehection group 
is the symmetric group ©„ of permutations of an n-element set. In unni, there appeared 
two determinant evaluations associated to the symmetric group, see Theorems 55 and 
56 in jlU9j . They concerned evaluations of determinants of the type 

det (^ (5.87) 

CTjTrGSn V / 

due to Varchenko, Zagier, and Thibon, respectively, in which stat is the statistic “num¬ 
ber of inversions, ” respectively “major index. ” We know that in many helds of math¬ 
ematics there exist certain diseases which are typical for that held. Algebraic combi¬ 
natorics is no exception. Here, I am not talking of the earlier mentioned “g-disease” 
(see Footnote |2ni although, due to the presence of q we might also count it as a case of 
g-disease), but of the disease which manifests itself by the question “And what about 
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the other types?” So let us ask this question, that is, are there theorems similar to 
the two theorems which we mentioned above for other reflection groups? 

So, hrst of all we need analogues of the statistics “number of inversions” and “major 
index” for other reflection groups. Indeed, these are available in the literature. The 
analogue of “number of inversions” is the so-called length of an element in a Coxeter 
group (see [HI] for the dehnition). As a matter of fact, a closed form evaluation of the 
determinant ( 071 ) . where ©„ is replaced by any hnite or affine reflection group, and 
where stat is the length, is known (and was already implicitly mentioned in |l()9j j. This 
result is due to Varchenko [TMl Theorem (1.1), where a{H) is specialised to q]. His 
result is actually much more general, as it is valid for real hyperplane arrangements in 
which each hyperplane is assigned a different weight. I will not state it here explicitly 
because I do not want to go through the dehnitions and notations which would be 
necessary for doing that. 

So, what about analogues of the “major index” for other reflection groups? These 
are also available, and there are in fact several of them. The hrst person to introduce a 
major index for rehection groups other than the symmetric groups was Reiner in USD]. 
He proposed a major index for the hyperoctahedral group Bn, which arose naturally 
in his study of P-partitions for signed posets. The elements of Bn are often called 
signed permutations, and they are all elements of the form 7ri7r2 ... 7r„, where Tj G 
{±1, ±2, ..., ±n}, i = 1,2,... ,n, and where |vri||7r2| ... |vr„| is a permutation in G„. To 
dehne their multiplicative structure, it is most convenient to view tt = 7ri7r2... vr„ as 
a linear operator on M” acting by permutation and sign changes of the co-ordinates. 
To be precise, the action is given by 7i{ei) = (sgn7ri)e|^.|, where is the Tth standard 
basis vector in M"’, i = 1,2,... ,n. The multiplication of two signed permutations is 
then simply the composition of the corresponding linear operators. 

The major index maj^ tt of an element 7i E Bn which Reiner dehned is, as in the 
symmetric group case, the sum of all positions of descents in vr. (There is a natural 
notion of “descent” for any Coxeter group.) Concretely, it is 


n—1 

maj^TT := x(7r„ <0) + '^i- xi^^i >b T^i+i), 

i=l 


where we impose the order 1 <b 2 <b ■ ■ ■ <b n <b —n • • • <b —2 —1 on our 

ground set {±1, ±2,..., ±n}, and where x{^) = 1 if .4, is true and x(w4,) = 0 otherwise. 


order to give a reader who is not acquainted with the language and theory of reflection groups 
an idea what this question is referring to, I mention that all finite reflection groups have been classified, 
each having been assigned a certain “type.” So, usually one proves something for the symmetric group 
&n, which, according to this classification, has type A„_i, and then somebody (which could be oneself) 
will ask the question “Can you also do this for the other types?”, meaning whether or not there exists 
an analogous result for the other finite reflection groups. 
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There is overwhelming computational evidence^^ that the “major-determinant” for 
Bn, i.e., the determinant (071) with stat = maj^ and with replaced by Bn, factors 
completely into cyclotomic polynomials. 


Conjecture 46. For any positive integer n, we have 

n n 

i=l i=2 

A different major index for Bn was proposed by Adin and Roichman in [S]. It arises 
there naturally in a combinatorial study of polynomial algebras which are diagonally 
invariant under Bn- (In fact, more generally, wreath products of the form Cml&n, where 
Cm is the cyclic group of order m, and their diagonal actions on polynomial algebras 
are studied in jS] . These groups are also sometimes called generalised reflection groups. 
In this context. Bn is the special case C 2 I &n-) If we write neg tt for the number of i for 
which TTi is negative, then the flag-major index fmaj of Adin and Roichman is dehned 
by 

fmaj vr := 2 majyi TT-|-negTT, (5.89) 

where maj^ is the “ordinary” major index due to MacMahon, 

n—1 

maj^TT ■=^i- x(vri > VTi+i). 

i=l 

If one now goes to the computer and calculates the determinant on the left-hand side of 
()5.88|1 with majg replaced by fmaj for n = 1,2, 3,4, 5 (see Footnote |27| for the precise 

^^The reader may wonder what this computational evidence could actually be. After all, we are 
talking about a determinant of a matrix of size 2"'n!. More concretely, for n = 1,2, 3,4, 5 these are 
matrices of size 2, 8, 48, 384, 3840, respectively. While Maple or Mathematica have no problem 
to compute these determinants for n = 1 and n = 2, it takes already considerable time to do the 
computation for n = 3, and it is, of course, completely hopeless to let them compute the one for n = 4, 
a determinant of a matrix of size 384 which has polynomial entries (cf. Footnote II811 . However, the 
results for n = 1,2,3 already “show” that the determinant will factor completely into factors of the 
form 1 — gh j = 1,2,..., 2n. One starts to expect the same to be true for higher n. To get a formula 
for n = 4, one would then apply the tricks explained in Footnote ^1 That is, one specialises q to 4, 
at which value the first 8 cyclotomic polynomials (in fact, even more) are clearly distinguishable by 
their prime factorisations, and one computes the determinant. The exponents of the various factors 
1 — g* can then be extracted from the exponents of the prime factors in the prime factorisation of 
the determinant with g = 4. Unfortunately, the data collected for n = 1,2, 3,4 do not suffice to 
come up with a guess, and, on the other hand. Maple and Mathematica will certainly be incapable to 
compute a determinant of a matrix of size 3840 (which, just to store it on the disk, occupies already 10 
megabytes ...). So then, what did I mean when I said that the conjecture is based on data including 
n = 5? This turned out to become a “test case” for LinBox, a C-I--I- template library for exact high- 
performance linear algebra Ea, which is freely available under http://linalg.org. To be honest, I 
was helped by Dave Saunders and Zhendong Wan (two of the developers) who applied LinBox to do 
rank and Smith normal form computations for the specialised matrix with respect to various prime 
powers (each of which taking several hours). The specific computational approach that worked here is 
quite recent (thus, it came just in time for our purpose), and is documented in m- The results of 
the computations made it possible to come up with a “sure” prediction for the exponents with which 
the various prime factors occur in the prime factorisation of the specialised determinant. As in the 
case n = 4, the exponents of the various factors 1 — g®, z = 1,2,..., 2n can then easily be extracted. 
(The guesses were subsequently also tested with special values of g other than g = 4.) 


det ( 

cr.TT^Bn \ 
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meaning of “calculating the determinant for n = 1,2, 3,4, 5”), then again the results 
factor completely into cyclotomic factors. Even more generally, it seems that one can 
treat the two parts on the right-hand side of (j5.89|l . that is “major index” and “number 
of negative letters,” separately. 


Conjecture 47. For any positive integer n, we have 


det 

(TjTrGBn 


^maj^{o-7r ^)j^neg{an l) 


i=l i=2 


(5.90) 


I should remark that Adin and Roichman have shown in [H] that the statistics fmaj 
is equidistributed with the statistics length on Bn- However, even in the case where we 
just look at the flag-major determinant (that is, the case where g = in ConiecturelTfll. 
this does not seem to help. (Neither length nor flag-major index satisfy a simple law 
with respect to multiplication of signed permutations.) In fact, from the data one sees 
that the flag-major determinants are different from the length determinants (that is, 
the determinants (071) . where 6„ is replaced by Bn and stat is flag-major, respectively 
length). 

Initially, I had my program wrong, and, instead of taking the (ordinary) major index 
maj^ of the signed permutation tt = 7ri7r2 ... 7r„ in (I5.89jl . I computed taking the major 
index of the absolute value of tt. This absolute value is obtained by forgetting all signs 
of the letters of tt, that is, writing |7r| for the absolute value of tt, |7r| = |7ri| |7r2| ... Ivr^,!. 
Curiously, it seems that also this “wrong” determinant factors nicely. (Again, the 
evidence for this conjecture is based on data which were obtained in the way described 
in Footnote 1771 j 


Conjecture 48. For any positive integer n, we have 


det^ _ ^2^2" in-n! _ ^i^2"n!(i-l)/i_ 


(7,7V£Br, 


i=2 


(5.91) 


Since, as I indicated earlier, Adin and Roichman actually dehne a flag-major index for 
wreath products Cml&n, a question that suggests itself is whether or not we can expect 
closed product formulae for the corresponding determinants. Clearly, since we are now 
dealing with determinants of the size rn^nl, computer computations will exhaust our 
computer’s resources even faster if m > 2. The calculations that I was able to do suggest 
strongly that there is indeed an extension of the statement in Conjecture 07] to the case 
of arbitrary m, if one uses the dehnition of major index and the “negative” statistics 
for Cm I &n as in jS]. (See 0 Section 3] for the dehnition of the major index. The sum 
on the right-hand side of jSl (3.1)] must be taken as the extension of the “negative” 
statistics neg to Cm I &n-) 


Problem 49. Find and prove the closed form evaluation of 

det {^ 

cr,-K£Cml&n. V / 


(5.92) 


where maj and neg are the extensions to Cm I &n of the statistics maj^ and neg in 
Conjecture ETl as described in the paragaph above. 


Together with Brenti, Adin and Roichman proposed another major statistics for 
signed permutations in They call it the negative major index., denoted nmaj, and 
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it is defined as the sum of the ordinary major index and the sum of the absolute values 
of the negative letters, that is, 

nmaj vr := maj^ vr + sneg tt, 

where snegvr := — < OjTTj. Also for this statistics, the corresponding de¬ 

terminant seems to factor nicely. In fact, it seems that one can again treat the two 
components of the dehnition of the statistics, that is, “major index” and “sum of nega¬ 
tive letters,” separately. (Once more, the evidence for this conjecture is based on data 
which were obtained in the way described in Footnote 1771 ) 


Conjecture 50. For any positive integer n, we have 


det 

cr,TrGBn 


^gmaj^(cr7r l)pSneg(o-7r 



i=l i=2 


(5.93) 


If one compares the (conjectured) result with the (conjectured) one for the “flag- 
major determinant” in Conjecture 071 then one notices the somewhat mind-boggling 
fact that one obtains the right-hand side of ()5.93|1 from the one of (I5.9()j) by simply 
replacing (in the factored form of the latter) 1 — by 1 — everything else, the 
exponents, the “g-part”, is identical. It is difficult to imagine an intrinsic explanation 
why this should be the case. 

Since Thibon’s proof of the evaluation of the determinant with stat being the 

(ordinary) major index for permutations (see jl()9[ Appendix C]) involved the descent 
algebra of the symmetric group, viewed in terms of non-commutative symmetric func¬ 
tions, one might speculate that to prove Conjectures ITHUTHI andlHniit may be necessary to 
work with Bn versions of descent algebras (which exist, see jl75j j and non-commutative 
symmetric functions (which also exist, see jS2])- 

Adriano Garsia points out that all the determinants in (15.8711 . (I5.88jl . (I5.9()|l - (l5.93jl 
are special instances of group determinants. (See the excellent survey article |12()j 
for information on group determinants.) The main theorem on group determinants, 
due to Frobenius, says that a general group determinant factorises into irreducible 
factors, each of which corresponding to an irreducible representation of the group, 
and the exact exponent to which the irreducible factor is raised is the degree of the 
corresponding irreducible representation. In view of this, Garsia poses the following 
problem, a solution of which would rehne the above conjectures. Problem 03 and the 
earlier mentioned results of Varchenko, Zagier, and Thibon. 


Problem 51. For each of the above special group determinants, determine the closed 
formula for the value of the irreducible factor corresponding to a fixed irreducible rep¬ 
resentation of the involved group {&n, Bn, Cm I &n, respectively). 

It seems that a solution to this problem has not even been worked out for the de¬ 
terminant which is the subject of the results of Varchenko and Zagier, that is, for the 
determinant with stat being the number of inversions. 

For further work on statistics for (generalised) reflection groups (thus providing fur¬ 
ther prospective candidates for forming interesting determinants), I refer the reader to 
[Iii3ii3i23i2nii2nii27iiniiis3isnii7niiiii. i must report that, somewhat disappoint¬ 
ingly, it seems that the various major indices proposed for the group Dn of even signed 
permutations (see iBiiziiinn]) apparently do not give rise to determinants in the same 
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way as above that have nice product formulae. This remark seems to also apply to 
determinants formed in an anlogous way by using the various statistics proposed for 
the alternating group in inHi. 

5.8. More poset and lattice determinants. Continuing the discussion of de¬ 
terminants which arise under the influence of the above-mentioned “reflection group 
disease,” we turn our attention to two miraculous determinants which were among the 
last things Rodica Simion was able to look at. Some of her considerations in this 
direction are reported in jlbhj . 

The first of the two is a determinant of a matrix the rows and columns of which are 
indexed by type B non-crossing partitions. This determinant is inspired by the eval¬ 
uation of an analogous one for ordinary non-crossing partitions (that is, in “reflection 
group language,” type A non-crossing partitions), due to Dahab j^O] and Tutte |18dj 
(see jlU9L Theorem 57, (3.69)]). Recall (see |177l Ch. 1 and 3] for more information) 
that a (set) partition of a set S' is a collection {Bi, B 2 , ..., B^] of pairwise disjoint non¬ 
empty subsets of S such that their union is equal to S. The subsets B^ are also called 
blocks of the partition. One partially orders partitions by refinement. With respect to 
this partial order, the partitions form a lattice. We write vr 7 (the A stands for the 
fact that, in “reflection group language”, we are looking at “type A partitions”) for 
the join of n and 7 in this lattice. Roughly speaking, the join of n and 7 is formed by 
considering altogether all the blocks of tt and 7 . Subsequently, whenever we hnd two 
blocks which have a non-empty intersection, we merge them into a bigger block, and 
we keep doing this until all the (merged) blocks are pairwise disjoint. 

If S' = { 1 , 2 ,... ,n}, we call a partition non-crossing if for any i < j < k < I the 
elements i and k are in the same block at the same time as the elements j and I are in 
the same block only if these two blocks are the same. (I refer the reader to jl 68 j for a 
survey on non-crossing partitions.) 

Reiner HSH introduced non-crossing partitions in type B. Partitions of type Bn are 
(ordinary) partitions of (1, 2,..., n, —1, —2,..., —n} with the property that whenever 
R is a block then so is —B := {—h : b G B}, and that there is at most one block B with 
B = —B. A block B with B = —B, if present, is called the zero-block of the partition. 
We denote the set of all type R„ partitions by fl^, the number of zero blocks of a 
partition vrby zbkvr, and we write nzbkvr for half of the number of the non-zero blocks. 
Type Bn non-crossing partitions are a subset of type R„ partitions. Imposing the order 
1<2<---<?7,<—1<—2<---< —n on our ground-set, the dehnition of type 
Bn non-crossing partitions is identical with the one for type A non-crossing partitions, 
that is, given this order on the ground-set, a R„ partition is called non-crossing if for 
any i < j < k < I the elements i and k are in the same block at the same time as the 
elements j and I are in the same block only if these two blocks are the same. We write 
NC^ for the set of all Bn non-crossing partitions. 

The determinant dehned by type R„ non-crossing partitions that Simion tried to 
evaluate was the one in (OUl below.The use of the “type A” join for two 

fact, instead of Va, the “ordinary” join, she used the join in the type Bn partition lattice 11^. 
However, the number of non-zero blocks will be the same regardless of whether we take the join of 
two type Bn non-crossing partitions with respect to “ordinary” join or with respect to “type i?„” join. 
This is in contrast to the numbers of zero blocks, which can differ largely. (To be more precise, one 
way to form the “type Bn’’ join is to first form the “ordinary” join, and then merge all zero blocks 
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type B non-crossing partitions in (j5.94ll may seem strange. However, this is certainly a 
well-defined operation. The result may neither be a type B partition nor a non-crossing 
one, it will just be an ordinary partition of the ground-set {1,2,..., n, —1, —2,..., —n}. 
We extend the notion of “zero block” and “non-zero block” to these objects in the 
obvious way. Simion observed that, as in the case of the analogous type An partition 
determinant due to Dahab and Tutte, it factors apparently completely into factors 
which are Chebyshev polynomials. Based on some additional numerical calculations,^® 
I propose the following conjecture. 


Conjecture 52. For any positive integer n, we have 


n 


det 

7r,7eNC® 


bk(7rVA7) 


)=n 

i=l 


/ t/3i-i(Vg/2) \(«-.) 


(5.94) 


where Um{x) := '^) (2^)™ m-th Chebyshev polynomial of the 

second kind. 


If proved, this would solve Problem 1 in |lb5j . It would also solve Problem 2 from 
[165j . because [/m-i(y^/2) is, up to multiplication by a power of g, equal to the product 
Y\j\rafA^)^ where the polynomials /j(g) are the ones of jl65j . A simple computation 
then shows that, when the right-hand side product of (jEMD is expressed in terms of 
the /;(g)’s, one obtains 


where 

This agrees with the data in jlfibj 

Footnote 1221) • 

into one big block.) The reason that I insist on using \Ja is that this is crucial for the more general 
Conjecture ISSl To tell the truth, the discovery of the latter conjecture is due to a programming error 
on my behalf (that is, originally I aimed to program the “type B” join, but it happened to be the 
“ordinary” join ...). 

^®Evidently, more than five years later, thanks to technical progress since then, one can go much 
farther when doing computer calculations. The evidence for Conjecture 1521 which I have is based 
on, similar to the conjectures and calculations on determinants for signed permutations in Subsec¬ 
tion o (see Footnote EH), the exact form of the determinants for n = 1,2, 3,4, which were already 
computed by Simion, and, essentially, the exact form of the determinants for n = 5 and 6. By “es¬ 
sentially” I mean, as earlier, that I computed the determinant for many special values of q, which 
then let me make a guess on the basis of comparison of the prime factors in the factorised results 
with the prime factors of the candidate factors, that is the irreducible factors of the Chebyshev 
polynomials. Finally, for guessing the general form of the exponents, the available data were not 
sufficient for Rate (see Footnote m. So I consulted the fabulous On-Line Encyclopedia of Integer Se¬ 
quences (http://www.research.att.com/~njas/sequences/Seis.html), originally created by Neil 
Sloane and Simon Plouffe nza, and since many years continuously further developed by Sloane and 
his team CZH An appropriate selection from the results turned up by the Encyclopedia then led to 
the exponents on the right-hand side of IS22I). 


n 

\{hk{qr-'\ (5.95) 

k=l 



(mod 3) 

and with the further ones I have computed (see 
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Even more seems to be true. The following conjecture predicts the evaluation of the 
more general determinant where we also keep track of the zero blocks. 

Conjecture 53. For any positive integer n, we have 

det (g-bkKv^7)^zbk(WA7)) = IT ( 2 T 2 i{^/ 2 ) + 2 - z) (5.96) 

where Tm{x) := m-th Chebyshev polynomial of 

the first kind. 


Again, the conjecture is supported by extensive numerical calculations. It is not too 
difficult to show, by using some identities for Chebyshev polynomials, that Conjec¬ 
ture ESI implies Conjecture ESI 

The other determinant which Simion looked at (cf. jl65[ Problem 9ff]), was the Bn 
analogue of a determinant of a matrix the rows and columns of which are indexed 
by non-crossing matchings, due to Lickorish [Ea, and evaluated by Ko and Smolinsky 
mn] and independently by Di Francesco [SHI (see [TIM Theorem 58]). As we may regard 
(ordinary) non-crossing matchings as partitions all the blocks of which consist of two 
elements, we define a Bn non-crossing matching to be a Bn non-crossing partition all the 
blocks of which consist of two elements. We shall be concerned with i? 2 n non-crossing 
matchings, which we denote by NCmatch(2n). With this notation, the following seems 
to be true. 

Conjecture 54. For any positive integer n, we have 

det tgnzbk(w^7)^^bk{w^7)') = IT (2T2i{q/2) + 2-z^\ . (5.97) 

7r,7€NCmatch(2ri.) \ / 

2=1 

The reader should notice the remarkable fact that, in the case that Conjectures IKHl and 
l^are true, the right-hand side of (ICTI) is, up to a power of z, equal to the right-hand 
side of ()5.96|1 with q replaced by q^ and z replaced by z^. An intrinsic explanation why 
this should be the case is not known. An analogous relation between the determinants 
of Tutte and of Lickorish, respectively, was observed, and proved, in BH. Also here, no 
intrinsic explanation is known. 

The reader is referred to |165j for further open problems related to the determinants 
in ConjecturesFinally, it may also be worthwhile to look at determinants dehned 
using Dn non-crossing partitions and non-crossing matchings, see and cni for two 
possible dehnitions of those. 

The reader may have wondered why in Con j ectures HM and 1531 we considered determi¬ 
nants dehned by type Bn non-crossing partitions, which form in fact a lattice, but used 
the extraneous type A join^° in the dehnition of the determinant, instead of the join 
which is intrinsic to the lattice of type Bn non-crossing partitions. In particular, what 
would happen if we would make the latter choice? As it turns out, for that situation 
there exists an elegant general theorem due to Lindstrom |125j . which I missed to state 
in |109j . I refer to CHI Ch. 3] for the explanation of the poset terminology used in the 
statement. 

^*^and not even the one in the type A non-crossing partition lattice! 






62 


C. KRATTENTHALER 


Theorem 55. Let L be a finite meet semilattice, R be a commutative ring, and f : 
L X L ^ R be an incidence function, that is, f{x, y) = 0 unless x Ay = x. Then 


detL{f{xAy,x)) = J] I y{x,y)f{x,y) 

y&L \x£L 


where y is the Mobius function of L. 


(5.98) 


□ 


Clearly, this does indeed answer our question, we just have to specialise f{x, y) = h{x) 
for X Ay = X, where h is some function from L to R. The fact that the above theorem 
talks about meets instead of joins is of course no problem because this is just a matter 
of convention. 

Having an answer in such a great generality, one is tempted to pose the problem of 
hnding a general theorem that would encompass the above-mentioned determinant eval¬ 
uations due to Tutte, Dahab, Ko and Smolinsky, Di Francesco, as well as Conjectures |S21 
and ESI This problem is essentially Problem 6 in |lb5j . 

Problem 56. LetL and L' be two lattices {semilattices?) with L' C L. Furthermore, let 
R be a commutative ring, and let f be a function from L to R. Under which conditions 
is there a compact formula for the determinant 


det {f{xALy)), 

x,yGL 


(5.99) 


where Al is the meet operation in L? 


By specialisation in Theorem 1551 one can derive numerous corollaries. For example, 
a very attractive one is the evaluation of the “GCD determinant” due to Smith na. 
(In fact. Smith’s result is a more general one for factor closed subsets of the positive 
integers.) 

Theorem 57. For any positive integer n, we have 


i=i 


where denotes the Euler totient function. 


□ 


An interesting generalisation of Theorem EHl to posets was given by Altini§ik, Sagan 
and Tuglu [Hj. Again, all undehned terminology can be found in mil ch. 3]. 

Theorem 58. Let P be a finite poset, R be a commutative ring, and f,g:PxP^R 
be two incidence functions, that is, f{x,y) = 0 unless x <y in P, the same being true 
for g. Then 


det 

x,yeP 


'^fiz,x)g{z,y) = Ylfix,x)g{x,x). 


\zeP 


xGP 


□ 


The reader is referred to Section 3 of [HI for the explanation why this theorem implies 
Lindstrom’s. 
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5.9. Determinants for compositions. Our next family of determinants consists of 
determinants of matrices the rows and columns of which are indexed by compositions. 
Recall that a composition of a non-negative integer n is a vector {ai,a 2 , ■ ■ ■ ,ak) of 
non-negative integers such that ai + a 2 + ■ ■ ■ + au = n, for some k. For a hxed 
fc, let C{n,k) denote the corresponding set of compositions of n. While working on 
a problem in global optimisation, Brunat and Montes HOI discovered the following 
surprising determinant evaluation. It allowed them to show how to explicitly express a 
multivariable polynomial as a difference of convex functions. In the statement, we use 
standard multi-index notation; if a = (ai, 0 : 2 , ■ ■ ■, Ufc) and (3 = (/?!, /d 2 , • • •, Pk) are two 
compositions, we let 


OL^ := 



Theorem 59. For any positive integers n and k, we have 

n—1 

det (af) = „(“L‘)+'=-i 

o:,/3€C(n,/c) 


(5.100) 

□ 


In recent joint work 123 EH, Brunat, Montes and the author showed that there is in 
fact a polynomial generalisation of this determinant evaluation. 

Theorem 60. Let x = (xi,a; 2 ,... ,Xk) he a vector of indeterminates, and let X be an 
indeterminate. Then, for any non-negative integers n and k, we have 

n 

det ((x + Xocf ) = (|x| + An)(D‘-) (5.101) 

i=l 

where x -|- Ao: is short for (xi -|- Aoi, X 2 + Aa 25 ■ ■ ■ ,Xk + Xak), and where |x| = xi + X 2 + 

-hXfc. 


As a matter of fact, there is actually a binomial variant which implies the above 
theorem. Extending our multi-index notation, let 

/aA _ /oiA / a2\ fak 
\(3) \(3j\(32)"\f3k 


Theorem 61. Let x = [xi,X 2 ,... ,Xk) he a vector of indeterminates, and let A be 
an indeterminate. Then, using the notation from Theorem inUl for any non-negative 
integers n and k, we have 


det 

a,l3(^C{n,k) 


X -|- Ack 

/3 


= Yi 


2 = 1 


-h (A - l)n -h i 


T,-|-fe —f—lA 
fe-1 ) 


(5.102) 


The above theorem is proved in [ 22 ] by identification of factors (see “Method 3” in 
Section ID). Theorem inn follows by extracting the highest homogeneous component in 
()5.101|1 . I report that, if one naively replaces “compositions” by “integer partitions” 
in the above considerations, then the arising determinants do not have nice product 
formulae. 


Another interesting determinant of a matrix with rows indexed by compositions ap¬ 
pears in the work of Bergeron, Reutenauer, Rosas and Zabrocki EH Theorem 4.8] on 
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Hopf algebras of non-commutative symmetric functions. It was used there to show that 
a certain set of generators of non-commutative symmetric functions were algebraically 
independent. To state their determinant evaluation, we need to introduce some nota¬ 
tion. Given a composition a = {ai,a 2 , ■ ■ ■ ,ak) of n with all summands positive, 
we let D{cil) = {ai, ai -|- 02 ,..., ai -f Q !2 + ■ ■ ■ + ttfc-i}- Furthermore, for two com¬ 
positions CK and f3 of n, we write oc U (3 for the (unique) composition 7 of n with 
D{‘y) = D{ct) U D{(3). Finally, a! is short for af. a 2 \ ■ ■ ■ ttfc!. 

Theorem 62. Let Comp(? 7 ,) denote the set of all compositions of n all summands of 
which are positive. Then 

e(-i) 

d^et ((aU/3)!)= JJ (5.103) 

Q;,pEComp(n) 

7 €Comp(n) i=l 

where i{'y) is the number of summands {components) of the composition 7 , and where 
Om denotes the number 0 /indecomposable permutations of m {cf. |178[ Ex. 5.13(b)]). 
These numbers can be computed recursively by ai = 1 and 

n —1 

a„ = n! — aj(n — 1)!, n > 1. 

i=l 

□ 

As explained in |2I] , one proves the theorem by factoring the matrix in (I5.1()3j) in the 
form CDC^, where C is the “incidence matrix” of “rehnement of compositions,” and 
where H is a diagonal matrix. Thus, in particular, the LU-factorisation of the matrix 
is determined. 

5.10. Two PARTITION DETERMINANTS. On the surface, integer partitions (see below 
for their dehnition) seem to be very closely related to compositions, as they can be 
considered as “compositions where the order of the summands is without importance.” 
However, experience shows that integer partitions are much more complex combina¬ 
torial objects than compositions. This may be the reason that the “composition de¬ 
terminants” from the preceding subsection do not seem to have analogues for integer 
partitions. Leaving aside this disappointment, here is a determinant of a matrix in 
which rows and columns are indexed by integer partitions. This determinant arose in 
work on linear forms of values of the Riemann zeta function evaluated at positive inte¬ 
gers, although the traces of it have now been completely erased in the hnal version of 
the article [TTIj . (The symmetric function calculus in Section 12 of the earlier version 
imi gives a vague idea where it may have come from.) 

Recall that the power symmetric function of degree d in xi,X 2 ,... ,Xk is given by 
xf X 2 + ■ ■ ■ + xf, and is denoted by Pd{xi, X 2 ,..., Xk). (See jl 21 [ Ch. 1 and 2 ], |13()l 
Ch. I] and jl781 Ch. 7] for in-depth expositions of the theory of symmetric functions.) 
Then, while working on dl, Rivoal and the author needed to evaluate the determinant 

det (pA(hi,h 2 , • • • ,hfc)), (5.104) 

A,/i€Part(n,/c) 

where Part(n, k) is the set of integer partitions of n with at most k parts, that is, the 
set of all possibilities to write n as a sum of non-negative integers, n = Ai-|-A 2 -I-- ■ - -l-Afc, 
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with Ai > A 2 > ■ ■ ■ > Afc > 0 (the non-zero Aj’s being called the parts of A), and where 
Pxixi,X2, ...,Xk)= Pxiixi,.. .,Xk)pxA^i, ...,Xk)-- ■Px^.ixi,.. .,Xk). 

Following the advice given in Section |2l we went to the computer and let it calculate 
the prime factorisations of the values of this determinant for small values of n and k. 
Indeed, the prime factors turned out be always small so that we were sure that a “nice” 
formula exists for the determinant. However, even more seemed to be true. Recall that, 
in order to facilitate a proof of a (still unknown) formula, it is (almost) always a good 
idea to try to introduce more parameters (see [TIM Sec. 2]). This is what we did. It 
led us consider the following determinant, 

det {px{pii + Xi, pi 2 + X 2 ,..., tik + Xk)), (5.105) 

A,/i€Part(n,fc) 


where Xi,X 2 ,...,Xfc are indeterminates. Here are some values of the determinant 


(nme 

for special values of n and k 

. For 

n 

= k = 

3 we obtain 




6 (Xi - X2 

+ 2 )(Xi-X 3 

+ 1)(X2 

-X3 

+ l)(Xi+X2 

+ X3 + 

3 )^ 

for n = 

4 and k = 3 we obtain 







8(Xi 

-X2 + l)(Xi 

-X2 + 2 )(Xi 

-X2 

+ 

3 )(A'i 

- X3 + 2) 








X 

to 

1 

CO 

+ 

+ 

to 

+ X3 + 4)" 

for n = 

5 and k = 3 we get 







8(Xi 

-X2 + 1)(X4 

-X2 + 2 )(X 4 

-X2 

+ 

3 )(Ai 

- X2 + 4 )(X 4 

1 

CO 

+ 

2) 


X (Xi - X 3 + 3 )(X 2 - X 3 + 1 )(X 2 - X 3 + 2)(Xi + X 2 + X 3 + 5)1\ 


for n = Q and k = 3 we get 

576 (Xi - X 2 + l)(Xi - X 2 + 2)2(Xi - X 2 + 3)2(Xi - X 2 + 4)(Xi - X 2 + 5) 

X (Xi - X 3 + l)(Xi - X 3 + 2)(Xi - X 3 + 3)(Xi - X 3 + 4)(X2 - X 3 + 1)2 

X (X 2 - X 3 + 2)2(Xi + X 2 + X 3 + 6 )'^ 

for n = 4 and A; = 4 we get 

192 (Xi - X 2 + l)(Xi - X 2 + 2)(Xi - X 2 + 3)(Xi - X 3 + 2)(X2 - X 3 + 1) 

X (Xi - X 4 + 1)(X2 - X 4 + 1)(X3 - X 4 + l)(Xi + X 2 + X 3 + X 4 + 4)^ 
while for n = 5 and fc = 4 we get 

48 (Xi - X 2 + l)(Xi - X 2 + 2)(Xi - X 2 + 3)(Xi - X 2 + 4)(Xi - X 3 + 2) 

X (Xi - X 3 + 3)(X2 - X 3 + 1)(X2 - X 3 + 2)(Xi - X 4 + 2) 

X (X 2 - X 4 + 1 )(X 3 - X 4 + l)(Xi + X 2 + X 3 + X 4 + 5 )^ 2 . 

It is “therefore” evident that there will be one factor which is a power of n + 
whereas the other factors will be of the form X* — Xj + Cij, again raised to some power. 
In fact, the determinant is easy to compute for k = 1 because, in that case, it reduces 
to a special case of the Vandermonde determinant evaluation. However, we were not 
able to come even close to computing ()5.1()5j) in general. (As I already indicated, hnally 
we managed to avoid the determinant evaluation in our work in HHO 
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To proceed further, we remark that evaluating (I5.1()5|) is equivalent to evaluating 
the same determinant, but with the power symmetric functions p\ replaced by the 
Schur functions Sx, because the transition matrix between these two bases of sym¬ 
metric functions is the character table of the symmetric group of the corresponding 
order (cf. |ld()[ Ch. I, Sec. 7]), the determinant of which is known (see Theorem IMI 
below; since in our determinants (I5.1()4j) and (Ib.lObj) the indices range over all par¬ 
titions of n with at most k parts, it is in fact the rehnement given in Theorem IHHl 
which we have to apply). The determinant with Schur functions has some advantages 
over the one with power symmetric functions since the former decomposes into a hner 
block structure. Alain Lascoux observed that, in fact, there is a generalisation of the 
Schur function determinant to Grafimannian Schubert polynomials, which contains an¬ 
other set of variables, Yi, Y 2 ,..., Yn+k-i- More precisely, given A = (Ai, A 2 ,..., Xk), let 
Ya(Xi, X 2 ,..., Xk,Yi,Y 2 ,...) denote the polynomial in the variables Xi, X 2 ,..., Xk 
and Yi,Y 2 , ..., dehned by (see inn Sections 1.4 and 9.7; the order of the Bk^ should 
be reversed in (9.7.2) and analogous places]) 


Yx{X,,X2,...,Xk;Yi,Y2,...) 


where the entries of the determinant are given by 


J2SmiXi,...,Xk;Yi,...,Yi)x^ 

m=0 




(5.106) 


The Grafimannian Schubert polynomial Ya reduces to the Schur function sx when all 
the variables T), i = 1, 2,..., are set equal to 0. Given these definitions, Alain Lascoux 
(private communication) established the following result. 


Theorem 63. Let Xi,X 2 ,..., Xk, Yi, Y 2 ,..., Yn^k-i be indeterminates. Then, 


det (Yx{pi + Xi, P 2 + X 2 ,..., pk + Xk, Yi, Y 2 ,...)) 

A,/iGPart(n,/c) 

n + 

a^G{n,k) \ 

where G{n, k) denotes the set of all Grafimannian permutations,^^ the descent of which 
{if existent) is at k, and which contain at most n — 1 inversions, and where n(?7,, k) is 


= n(n.fc) n 




a{i) 


(5.107) 


2=1 


2=1 


^^A permutation a in ©co (the set of all permutations of the natural numbers N which fix all but 
a finite number of elements of N) is called Grafimannian if a{i) < a{i + 1) for all i except possibly for 
one, the latter being called the descent of a (see [nni p. 13]). An inversion of cr is a pair {i,j) such 
that i < j but a{i) > (j{j). We remark that the number of Grafimannian permutations with descent 
(if existent) at k and at most n — 1 inversions is equal to the number of partitions of at most n — 1 
with at most k parts (including the empty partition). More concretely, if we denote this number by 
gn,k, the generating function of the numbers gn,k is given by 


gn,kT 

n—1 


n 


1 - 9 1 - 9* 
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given recursively by 

k 

n(n, k) = n(n, k — l)n(n — k, k) n n(v + /i,-Xfc), (5.108) 

a,rt^ {n,k) 

Part’'"(n, k) denoting the set of partitions ofn into exactly k (positive) parts, with initial 
conditions n(n, k) = 1 if k = 1 or n < 1. Explicitly, 

k-\w\B 

Tl(n,k)= JJ Yl Yl + hj - ^k-\w\B), (5.109) 

w fi£Pa,rt(n—\w\ji^k-\-m-vw,k—\w\B) j=l 
\w\B<k-l 

where the product over w runs over all finite-length words w with letters from {A, B}, 
including the empty word. The notation \w\b means the number of occurrences of B 
in w, with the analogous meaning for \w\a- The quantity invta denotes the number of 
inversions of w = WiW 2 ■ ■ ■, which is the number of pairs of letters (wi, Wj), i < j, such 
that Wi = B and Wj = A. 

We obtain the evaluation of the determinant (I5.104|l if we set X, = Y) = 0 for all i 
in the above theorem and multiply by (I5.114|l . Likewise, we obtain the evaluation of 
the determinant (j5.105|l if we set X* = 0 for all i in the above theorem and multiply the 
result by (I5.114j) . 

For example, here is the determinant (I5.107j) for n = 3 and k = 2, 

(Xi - X 2 + 2)(Xi + X 2 - Xi - X 2 + 3)(Xi + X 2 - Xi - X 3 + 3) 

X (Xi + X 2 - X 2 - X 3 + 3)(Xi + X 2 - Xi - X 4 + 3), 
and the following is the one for n = k = 3, 

(Xi - X 2 + 2)(Xi - X 3 + 1)(X2 - X 3 + 1) 

X (Xi + X 2 + X 3 - Xi - X 2 - X 3 + 3)(Xi + X 2 + X 3 - Xi - X 2 - X 4 + 3) 

X (Xi + X 2 + X 3 - Xi - X 3 - X 4 + 3)(Xi + X 2 + X 3 - Xi - X 2 - X 5 + 3). 

In the sequel, I sketch Lascoux’s proof of Theorem 1^21 The hrst step consists in 
applying the Monk formula for double Schubert polynomials in the case of Grafimannian 
Schubert polynomials (see EEl), 

/ k k \ 

- E = ’^^+(1,0,0,,,,) + E''<‘' 

\ 2=1 2=1 / fL 

where Ya is short for Ya(Xi, X 2 ,..., X*,; Xi, X 2 ,...), and where the sum on the right- 
hand side is over all partitions fv of the same size as A-|- (1, 0, 0,...) but lexicographically 
smaller. Clearly, by using this identity, appropriate row operations in the determinant 
()5.1()7|1 show that 

k k 

n-h ^Xj - (5.110) 

2=1 2=1 

is one of its factors, the relation between a and A being a(i) = Xk-i+i+i, i = 1,2,..., k. 
In particular, a can be extended (in a unique way) to a Grafimannian permutation. 
Moreover, doing these row operations, and taking out the factors of the form (I5.11()|l . 
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we collect on the one hand the product in (I5.1()7|l . and we may on the other hand 
reduce the determinant (j5.1()7|l to a determinant of the same form, in which, however, 
the partitions A = (Ai, A 2 ,..., Xk) run over all partitions of size at most n with the 
additional property that Ai = A 2 (instead of over all partitions from Part(n, k)). 

As it turns out, the determinant thus obtained is independent of the variables hi, ^ 2 , 
..., Yn+k-i- Indeed, if we expand each Schubert polynomial Ya(pi+Xi, /i 2 +-A 2 ,..., pfc+ 
Xk] hi, Y 2 ,...) in the determinant as a linear combination of Schur functions in Xi, X 2 , 
..., Xk with coefficients being polynomials in the hi, h 2 ,..., Yn+k-i, then, by also using 
that 

k 

+ Xi, p,2 + -A2,..., /ifc + Xk) = n Y Xi 

i=l 

is independent of /i, it is not difficult to see that one can eliminate all the h^’s by 
appropriate row operations. 

To summarise the current state of the discussion: we have already explained the oc¬ 
currence of the product on the right-hand side of (I5.1()7|) as a factor of the determinant. 
Moreover, the remaining factor is given by a determinant of the same form as in (I5.1()7|l . 


det (Ya(pi -|- Xi, p 2 + X 2 ,..., -|- Xk] hi, h 2 ,...)), (5.111) 

where, as before, p runs over all partitions in Part(n, k), but where A runs over partitions 
(Ai, A 2 ,..., Xk) of size at most n, with the additional restriction that Ai = A 2 . The hnal 
observation was that this latter determinant is in fact independent of the h^’s. This 
allows us to specify Yi arbitrarily, say hi = Xk- Now another property of double 
Schubert polynomials, namely that (this follows from the dehnition of double Schubert 
polynomials by means of divided differences, see | 121 [ (10.2.3)], and standard properties 
of divided differences) 


k 


Yp+(i,i,...,i)(Xi,X2,...,Xfc;hl,h'2,. 


Y,(Xi, X 2 , ...,Xk]Y2,...) 1[{X, - Yi) 


i=i 

(5.112) 

comes in handy. (In the vector (1,1,..., 1) on the left-hand side of ()5.112|1 there are k 
occurrences of 1.) Namely, if Yi = Xk, the matrix in ()5.111|1 (of which the determinant 
is taken), M{n, k] Yi, Y 2 ,...) say, decomposes in block form. If A is a partition with 
Afc > 0 and /i is a partition with = 0, then, because of (15.11211 . the corresponding 
entry in (I5.111|l vanishes. Furthermore, in the block where A and /i are partitions with 
Afc = /Ufc = 0, because of the dehnition (j5.1()0|l of the quantities Sm{-) the corresponding 
entry reduces to 


Ya(/^i + Xi, /i2 + X 2 , ..., p.fc-1 -|- Xk-i, Xk] Xk, Y 2 ,...) 

= + ^1) h2 + X 2 , . . . , /ifc_i -|- Xk-i] Y 2 , . . .). 

In other words, this block is identical with M{n, k — 1; Y 2 ,...). Finally, in the block 
indexed by partitions A and /i with A^ > 0 and /ifc > 0, we may use (I5.112|) to factor 
+ fij — Xk) out of the column indexed by p, for all such p. What remains 
is identical with M{n — k, k]Y 2 ,...). Taking determinants, we obtain the recurrence 
()5.1()8|1 . (Here we use again that the determinants in ()5.111|1 . that is, in particular, the 


























ADVANCED DETERMINANT CALCULUS: A COMPLEMENT 


69 


determinants of M(n, k; Yi, ...), M(n, fc — 1; 1^2, • • •): of M{n — k, k; Y 2 , ■ ■ ■), are all 

independent of the Yi^s.) The explicit form (I5.1()9|l for n(n, k) can be easily derived by 
induction on n and k. 

A few paragraphs above, we mentioned in passing another interesting determinant 
of a matrix the rows and columns of which are indexed by (integer) partitions: the 
determinant of the character table of the symmetric group 6„ (cf. (HHl Cor. 6.5]). Since 
this is a classical and beautiful determinant evaluation which I missed to state in uni, 
I present it now in the theorem below. There, the notation A h n stands for “A is a 
partition of n." For all undehned notation, I refer the reader to standard texts on the 
representation theory of symmetric groups, as for example jHi HOI I159j . 

Theorem 64. For partitions A and p ofn, letx^ip) denote the value of the irreducible 
character evaluated at a permutation of cycle type p. Then 


A^hn = n n (5.113) 

fL hn 2>1 

In words: the determinant of the character table of the symmetric group ©„ is equal to 
the products of all the parts of all the partitions of n. □ 


A rehnement of this statement, where we restrict to partitions of n with at most k 
parts, is the following. 


Theorem 65. With the notation of Theorem for all positive integers n and k, 
n > k, we have 


det 

A,p€Part(n,Ai) 


(x^(p)) 


pGPart(n,/i:) ^>1 


where mi{p) is the number of times i occurs as a part in the partition p. 


(5.114) 

□ 


This determinant evaluation follows from the decomposition of the full character 
table of ©,1 in the form 

= (5.115) 

where L = {Lx,fj,)\^fj,\-n is the transition matrix from power symmetric functions to 
monomial symmetric functions, and where K = (A'A,/i)A,Aii-n is the Kostka matrix, the 
transition matrix from Schur functions to monomial symmetric functions (see nsni 
Ch. I, (6.12)]). For, if we order the partitions of n so that the partitions in Part(n, k) 
come before the partitions in Part(n, k + 1), k = 1, 2,..., n — 1, and within Part(n, k) 
lexicographically, then, with respect to this order, the matrix L is lower triangular and 
the matrix K, and hence also K~^, is upper triangular. Furthermore, the matrix K, 
and hence also K~^, is even block diagonal, the blocks along the diagonal being the 
ones which are formed by the rows and columns indexed by the partitions in Part(n, k), 
k = 1,2,... ,n. These facts together imply that the decomposition (15.11511 restricts to 
the submatrices indexed by partitions in Part(n, k). 


(x^(p)) 


A,p€Part(n,fc) 


(-^A,p) A,pGPart(n,/i;) (Kx,) 


-1 

A,p.€Part(n,fc) * 


If we now take determinants on both sides, then, in view of = 1 and of 


i>l 
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for all yU, the theorem follows. 

Further examples of nice determinant evaluations of tables of characters of represen¬ 
tations of symmetric groups and their double covers can be found in 123 im. Deter¬ 
minants of tables of characters of the alternating group can be found in j24j . 

5.11. Elliptic determinant evaluations. In special functions theory there is cur¬ 
rently a disease rapidly spreading, generalising the earlier mentioned g-disease (see 
Footnote I2DI)- If could be called the “elliptic disease.” Recall that, during the g-disease, 
we replaced every positive integer nby l-|-g-|-g^-|-----|- g""^ = (1 — g"')/(l — g), 
and, more generally, shifted factorials a(a -|- 1) ■ ■ ■ (a -|- — 1) by g-shifted factorials 

(1 — «)(! — qa) ■ ■ ■ (1 — ag^”^). (Here, a takes the role of g“, and one drops the powers 
of 1 — g in order to ease notation.) Doing this with some “ordinary” identity, we ar¬ 
rived (hopefully) at its q-analogue. Now, once infected by the elliptic disease, we would 
replace every occurrence of a term 1 — x (and, looking at the dehnition of g-shifted 
factorials, we can see that there will be many) by its elliptic analogue 6{x]p): 

OO 

9{x) = 6{x]p) = ]^(1 —pf’x){l —pP^^/x). 

3=0 

Here, p is a complex number with \p\ < 1, which will be hxed throughout. Up to a 
trivial factor, equals the Jacobi theta function 6 i{x\t) (cf. jl89j h Clearly, 

9{x) reduces to 1 — t if p = 0. 

At hrst sight, one will be sceptical if this is a fruitful thing to do. After all, for 
working with the functions 9{x), the only identities which are available are the (trivial) 
inversion formula 

9(l/x) =--9(x), (5.116) 

X 

the (trivial) quasi-periodicity 

9{px) = ——9{x), (5.117) 

X 

and Riemann’s (highly non-trivial) addition formula (cf. jl89[ p. 451, Example 5]) 

u 

9ixy) 9(x/y) 9iuv) 9iu/v) — 9(xv) 9(x/v) 9(uy) 9(u/y) = — 9(yv) 9iy/v) 9(xu) 9ix/u). 

y 

(5.118) 

Nevertheless, it has turned out recently that a surprising number of identities from the 
“ordinary” and from the “g-world” can be lifted to the elliptic level. This is particularly 
true for series of hypergeometric nature. We refer the reader to Chapter 11 of [72] for 
an account of the current state of the art in the theory of, as they are called now, elliptic 
hypergeometric series. 

On the following pages, I give elliptic determinant evaluations a rather extensive 
coverage because, hrst of all, they were non-existent in |l()9j (with the exception of 
the mention of the papers by Milne), and, second, because I believe that the 

“elliptic research” is a research direction that will further prosper in the next future and 
will have numerous applications in many helds, also outside of just special functions 
theory and number theory. I further believe that the determinant evaluations presented 
in this subsection will turn out to be as fundamental as the determinant evaluations in 
Sections 2.1 and 2.2 in [innj. For some of them this belief is already a fact. For example. 
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determinant evaluations involving elliptic functions have come into the picture in the 
theory of multiple elliptic hyper geometric series, see PHl 11471115211155L11551117511187j . 
They have also an important role in the study of Ruijsenaars operators and related 
integrable systems IHH [mi. Furthermore, they have recently found applications in 
number theory to the problem of counting the number of representations of an integer 
as a sum of triangular numbers |154j . 

Probably the hrst elliptic determinant evaluation is due to Frobenius |nni (12)]. This 
identity has found applications to Ruijsenaars operators usi, to multidimensional el¬ 
liptic hypergeometric series and integrals ra. mzi and to number theory HMj. For a 
generalisation to higher genus Riemann surfaces, see jnni Cor. 2.19]. Amdeberhan [7] 
observed that it can be easily proved using the condensation method (see “Method 2” 
in Section ID). 

Theorem 66. Let Xi,X 2 , ■ ■ ■, Xn, Ui, 02 ,..., a„ and t be indeterminates. Then there 
holds 


/ ejta.x,) \ 

i^7<n ye{t) e{ajXi)) 


OjXj 9{ai/aj) 9{xi/xj) 

9(tai • • • QjiXi • • • Xjfj l<i<j<n 


9{t) 


n 

Li=i 


(5.119) 


□ 


For p = 0 and t —>■ 00 , this determinant identity reduces to Cauchy’s evaluation (ESD 
of the double alternant, and, thus, may be regarded as an “elliptic analogue” of the 
latter. 

Okada ma Theorem 1.1] has recently found an elliptic extension of Schur’s evalu¬ 
ation Q of a Cauchy-type Pfaffian. His proof works by the Pfaffian version of the 
condensation method. 


Theorem 67. Let Xi,X 2 , ■ ■ ■ ,Xn, t and w be indeterminates. Then there holds 

/ 9{xj/xi) 9{tXiXj) 9{wXiXj) \ 
l<i,j<2n y 9{XiXj) 9{t) 9{w) ) 

_9{tXi---X2n)9{wXi---X2n) TT 9[XjlXi) /..on', 

0(t) 0(u;) J^i^^^,^9ix^x,y ^ ^ 

□ 

The next group of determinant evaluations is from |1561 Sec. 3]. As the Vandermonde 
determinant evaluation, or the other Weyl denominator formulae (cf. |1()9[ Lemma 2]), 
are fundamental polynomial determinant evaluations, the evaluations in Lemma (HHl be¬ 
low are equally fundamental in the elliptic domain as they can be considered as the 
elliptic analogues of the former. Indeed, Rosengren and Schlosser show that they imply 
the Macdonald identities associated to affine root systems II2HI. which are the affine 
analogues of the Weyl denominator formulae. In particular, in this way they obtain 
new proofs of the Macdonald identities. 

In order to conveniently formulate Rosengren and Schlosser’s determinant evalua¬ 
tions, we shall adopt the following terminology from jl56j . For 0 < |p| < 1 and t y 0, 
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an An-i theta function f of norm t is a holomorphic function for x 7^ 0 such that 

/(px) = ^^/(x). (5.121) 

Moreover, if R denotes either of the root systems Bn, B'f, Cn, Cn, BCn or Dn (see 
Footnote |2ni and for information on root systems), we call f an R theta function if 



/(1/x) = — fix), 

X 

R = Bn, 

f(P^)- pn^2nf(^)^ 

/(1/x) = -fix). 

R = Bl, 

/(P^) - pn+1^2n+2 

/(1/x) = -fix). 

R = Cn, 


/(1/x) = --fix), 

X 

R = c:, 

f(P^) - pn^2n+l 

/(1/x) = --fix), 

X 

R = BCn, 

fiP^) - pn-1^2n-2 /(^)’ 

/(1/a^) = f{x). 

R = Dn. 

Given this dehnition, Rosengren and Schlosser jl56[ Lemma 3.2] show that a function 
/ is an An-i theta function of norm t if and only if there exist constants C, bi,... ,bn 
such that bi- ■ - bn = t and 

/(x) = C'6'(6ix) • ■■9{bnx), 


and for the other six cases, they show that / is an i? theta function if and only if there 
exist constants C, bi,..., bn-i such that 

fix) = C 9ix) 9ibix) 9ibi/x) ■ ■ ■ 9ib, 

^_ix) 9ibn-i/x), 

R = Bn, 

f{x) = C x~^9{x‘^-,p'^) 9{bix) 9{bi/x 

) ■ ■■9ibn-ix) 9ibn-i/x), 

R = Bl, 

f{x) = Cx~^9{x‘^) 9{bix) 9{bi/x) ■ ■ 

■9ibn-ix) 9ibn-i/x), 

R = Cn, 

/(x) = C9{x-,p2) 9{bix) 9{bi/x) ■ ■ ■ 

9ibn-\x) 9ibn-i/x), 

R = c:. 


/(x) = C9{x) 9{px‘^;p^) 9{bix) 6'(6i/x) ■ ■■9{bn-ix) 9{bn-i/x), R = BCn, 
f(x) = C 9{bix) 9{bi/x) ■ ■ ■ 9{bn-ix) 9{bn-i/x), R = Dn, 


where 9{x) = 9{x;p). 

If one puts p = 0, then an An-i theta function of norm t becomes a polynomial of 
degree n such that the reciprocal of the product of its roots is equal to t. Similarly, if 
one puts p = 0, then a Dn theta function becomes a polynomial in (x + l/x) of degree n. 
This is the specialisation of some of the following results which is relevant for obtaining 
the earlier Lemmas cnra 

The elliptic extension of the Weyl denominator formulae is the following formula. 
(See jl5fi[ Prop. 3.4].) 


Lemma 68. Let fi,---,fn be An-i theta functions of norm t. Then, 

det {fj {xi)) = C WAn-i (x), 


(5.122) 
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for some constant C, where 

Wa„_:,{x) = 9{tXi- ■ -Xn) JJ XjO{Xi/Xj). 

Moreover, if R denotes either Bn, Bf, Cn, Cf, BCn or Dn and /i,..., are R theta 
functions, we have 


, det {fj (a:i)) = C Wr{x), 

1<2 j<n 


for some constant C, where 


n Xi^9{xiXj)9{xi/xj), 

2=1 

n 

= Y[x-^9{xI,p^) JJ xf^9{xiXj)9{xi/xj), 

2=1 

n 

= Y[x-^9{x^i) JJ xf^9{xiXj)9{xi/xj), 

2=1 

n 

^c'v(t) = Jj6»(a:i;p5) JJ x~^9{xiXj) 9{xi/xj), 

2=1 

n 

Wbc„{x) = Y[9{xi)9{pxl,p^) JJ x-^9{xiXj)9{xi/xj), 
2=1 

= JJ xf^9{xiXj)9{xi/xj). 

l< 2 <ji'<n 


(5.123) 


□ 


Rosengren and Schlosser show in |15(il Prop. 6.1] that the famous Macdonald identi¬ 
ties for affine root systems jl28j are equivalent to special cases of this lemma. We state 
the corresponding results below. 


Theorem 69. The following determinant evaluations hold: 

{p-,p) 


det (x{-^9{{-l)^-^jf-Hx^;p^)) = , 


n 

oo 


Wa^_A^), 


l< 2 j<n 


Wb (^) 


det {xl ^xf^^p^'^) 

l<2j<n 




2(p";p")oo(p;p)^ ^ ^ . 

- Wsv(a;), 


(p2n.p; 


!2n\n 
oo 
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det 1 X' 

l< 2 ,jf<n 


j-n-l 


9{—p>x, 


272+2. ^272+2^ 


X, 


72+1-j 




{p-p) 


n 

oo 


^p2n+2. p2n+2' 


W'C.W. 


det^ (x^ "'9{—pP ^9{—pP 




(p2;p5)oo(p;p)^ .. 

■ W^cv(x), 


(p272.p: 


2n\n 
oo 


det (x^ '^9{—p’xf^^^]p^'^^^)—x^^^ ^9{—pPx^‘^^ 

l<i,j<n ' ' 


{p;p)[ 


^^271+1- 2^271-1-1 


WbcJx), 


,det (x’i "Oi-p’ +■” “;p’” ^)+i” +(-+ +■ “;p^” ^)) 

l<2j<n 


4 (p;p) 


n 

oo 


^p2n—2- p2n—2'^ 


—2 \n 
oo 


W^D„(x), n>2. 


□ 


Historically, aside from Frobenius’ elliptic Cauchy identity (Id.llbj) . the subject of 
elliptic determinant evaluations begins with Warnaar’s remarkable paper |187j . While 
the main subject of this paper is elliptic hypergeometric series, some elliptic determinant 
evaluations turn out to be crucial for the proofs of the results. Lemma 5.3 from HHH 
extends one of the basic determinant lemmas listed in uni, namely uni Lemma 5], 
to the elliptic world, to which it reduces in the case p = 0. We present this important 
elliptic determinant evaluation in the theorem below. 

Theorem 70. Let Xi,X 2 , ■ ■ ■, x„, Oi, 02 , ..., be indeterminates. For each j = 1,... ,n, 
let Pj{x) be a Dj theta function. Then there holds 


det \ Pj{xi) TT 9(akXi) 9(ak/xi 

1<2,7<72 \ 

- - ' k=j + l 


]^Pi(ai) QjXp 9{xjXi)9{xj/xi). (5.124) 

2=1 

□ 


Warnaar used this identity to obtain a summation formula for a multidimensional 
elliptic hypergeometric series. Further related applications may be found in |152l 11531 
The relevant special case of the above theorem is the following (see HHIl 
Cor. 5.4]). It is the elliptic generalisation of jlOOL Theorem 28]. In the statement, we 
use the notation 

{a-,q,p)m = 9{a-,p) 9{aq-,p) ■ ■ ■ 9{aq'^~^-,p), (5.125) 

which extends the notation for g-shifted factorials to the elliptic world. 
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Theorem 71. Let Xi, X2 ,..., A, B and C he indeterminates. Then, for any 
non-negative integer n, there holds 


l<^j<n\{BXi,q,p)n-j {BC/Xi,q,p)n-j) 



n X,0(X.IX,)0(CIX.X,)Y[ 

l<i<JI<n 2=1 


{B/A-, q, p)i_i {ABCq'^^-'^^] q, p)i-i 

{BXi-,q,p)n-i{BC/Xi-,q,p)n-i 

(5.126) 


□ 


Theorem 29 from [Tunj . which is slightly more general than [TIM Theorem 28], can also 
be extended to an elliptic theorem by suitably specialising the variables in Theorem 1701 

Theorem 72. Let Xi, X 2 ,..., X^, hi, Y 2 ,..., Y^, A and B be indeterminates. Then, 
for any non-negative integer n, there holds 

. t ( (X,Yf,q,p),{AY,/Xuq,p)A 
V (BA',; q,p), {AB/Xp q.p), ) 

= q^^^){ABp J] e{X,X,/A)e{X^/Xf) 

l<2<j<n 

^ A {ABYiq^-^-q,p)i_i {Yj/Bq^-^-q,p)i_i ^ 

\}^Xi-\BXp,q,p)^_^{AB/Xp,q,p)^.^' 

□ 


Another, very elegant, special case of Theorem 1701 is the following elliptic Cauchy- 
type determinant evaluation. It was used by Rains dm Sec. 3] in the course of deriving 
a. BCn BCjn integral transformation. 


Lemma 73. Let Xi, X 2 , ■ ■ ■, Xn and Oi, 02 ,..., a„ be indeterminates. Then there holds 


det 


6{ajXi) 9{aj/xi 


UjXj ^ 9{xjXi) 9{xj/xi) 9{aiaj) 9{ai/aj) 

n 

9{ajXi)9{aj/xi) 

ij=l 


(5.128) 


□ 

The remaining determinant evaluations in the current subsection, with the exception 
of the very last one, are all due to Rosengren and Schlosser |156j . The hrst one is a 
further (however non-obvious) consequence of Theorem 1701 fsee |156[ Cor. 4.3]). Two 
related determinant evaluations, restricted to the polynomial case, were applied in j162] 
and |164j to obtain multidimensional matrix inversions that played a major role in the 
derivation of new summation formulae for multidimensional basic hypergeometric series. 
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Theorem 74. Let Xi,X 2 , ■ ■ ■, Xn, ai, 02 ,, dn+i, b be indeterminates. For each 
j = 1,..., n + 1, let Pj{x) be a Dj theta function. Then there holds 


n +1 


Pn+iib) det \ Pj{xi) TT {0{akXi)6{ak/xi)) 

- - ' k=j+l 

n +1 


-fn+1 (^2 j 
Pn+l{b) 


Pj{b) J] {e{a,b)e{a,/b)) 


k=j+l 


n +1 


= Y\PMi) n 0(xjXi)9(xj/xi), (5.129) 

2=1 l<2<_7<n+l 

where Xn+i = b. □ 

The next determinant evaluation is Theorem 4.4 from |156j . It generalises another 
basic determinant lemma listed in uni, namely Lemma 6 from jl()9j . to the elliptic 
case. It looks as if it is a limit case of Warnaar’s in Theorem EOl However, limits are 
very problematic in the elliptic world, and therefore it does not seem that Theorem 1701 
implies the theorem below. For a generalisation in a different direction than Theorem 1701 
see jl82[ Appendix B] (cf. also [15(4 Remark 4.6]). 

Theorem 75. Let xi,X 2 , ■ ■ ■ ,Xn, ai, 02 ,..., and t be indeterminates. For each 
j = 1,.. . ,n, let Pj{x) be an Aj-i theta function of norm tai - ■ ■ aj. Then there holds 


^det^ Pj{xi) 0{akXi) 

' k=j+i 


6{tai ■ --anXi ■ - -Xn) 

e{t) 


Y\_Pi{,l/ai) ajXj9{xi/xj). (5.130) 

2=1 

□ 


As is shown in [1561 Cor. 4.8], this identity implies the following determinant evalu¬ 
ation. 

Theorem 76. Let xi,X 2 , ■■■ ,Xn, oi, 02 ,..., a^+i and b be indeterminates. For each 
j = 1,... ,n -f 1, let Pj{x) be an Aj_i theta function of norm tai ■ ■ - aj. Then there 
holds 


n +1 


Pn+i{b) det Pj{xi) TT 9{akXi) - 

l<i,j<n ' 


k=j+l 


Pfi+l (^ 2 ) 
Pn+lib) 


n+1 


Piib) n oM 


k=j+l 


9(tbai ■ --an+iXi ■■■Xn) 
9{t) 


n +1 


where Xn+i = b. 


]^Pi(l/ai) ajXj9{xi/xj), (5.131) 

2=1 l<2<_7<n+l 

□ 


By combining Lemma [701 and Theorem [701 a determinant evaluation similar to the 
one in Theorem [71 but different, is obtained in [1561 Theorem 4.9]. 
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Theorem 77 . Let Xi,X2, ■ ■ ■, Xn, Oi, 02,..., and Ci,..., Cn+2 be indeterminates. For 
each j = 1, ..., n, let Pj{x) be an Aj_i theta function of norm (ci • • • c„+2aj+i • • • 

Then there holds 


det 


n +2 


X, 


n 9 {akXi) 

k=j+l 

n-\-2 n 

-xf^^Pj{x-^)WO{cux-^) JJ 0{akXl 

k=j+l 


k=l 




k=l 


ai-■ ■ a„ 


9(^Ci • • • C^+2®1 ■ ■ ■ 


Y[Pi{l/ai) 


X J] 0 {ciCj)l[e{x^,) J] OjX^ ^ 9 {xiXj) 9 {xi/xj). 

l<i<j<n-\-2 i=l 


(5.132) 

□ 


The last elliptic determinant evaluation which I present here is a surprising elliptic 
extension of a determinant evaluation due to Andrews and Stanton m Theorem 8] 
(see |1091 Theorem 42]) due to Warnaar |1871 Theorem 4.17]. It is surprising because in 
the former there appear g-shifted factorials and g^-shifted factorials at the same time, 
but nevertheless there exists an elliptic analogue, and to obtain it one only has to add 
the p everywhere in the shifted factorials to convert them to elliptic ones. 

Theorem 78. Let x and y be indeterminates. Then, for any non-negative integer n, 
there holds 

f {(l/yxq"-,q‘^,p)i-j (l/x^g^+^b g^p)^-A 

0<io<n-i V (g; g, p)2i+i-j ( 1 / yxq^^; q, p)i-j (iz/xgi+^h g, p)i_j ) 

^ (T^g^^+^ g, p)i (V+Vi/; g^ p)i {yxq^+^; g^ p)i 

H (a;2g2i+2. q2^p^. q2^ {yxq^+^; g, p)i [xq^+^jy, q,p)i' 

□ 


In closing this hnal subsection, I remind the reader that, as was already said before, 
many Hankel determinant evaluations involving elliptic functions can be found in 11331 
and |mj . 
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